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Measurement  of  the  spatial  decay  of  spectral  amplitudes  of  the  higher-mode  seismic  surface  wave  train  L( 
indicates  that  in  the  Great  Basin,  the  apparent  seismic  quality  factor,  Q,  is  a  function  of  frequency.  Analysis  of 
recordings  from  a  NTS  nuclear  explosion  made  along  a  300  km  profile  through  northwest  Nevada  yield  the 
function  2(/>206/068  over  the  band  0.3</<10.0  Hz.  Similar  analysis  using  numerous  recordings  of  carth- 
quakcj  along  paths  which  provide  a  good  average  sampling  of  the  Great  Basin  crust  gives 
0(/)=214(±15)/054±009,  for  0.3</<5.0  Hz.  If  the  crustal  sampling  by  Lt  energy  from  nuclear  explosions  (sur¬ 
face  sources)  is  primarily  in  the  shallow  crust,  as  has  been  suggested,  then  these  results  indicate  a  greater  fre¬ 
quency  dependence  of  apparent  Q  there  than  at  depth.  To  better  estimate  the  errors  in  Q  at  each  frequency,  we 
have  used  a  moving  window  spectral  analysis  technique  described  by  Butler  et  al  (1987)  to  analysis  the  nuclear 
explosion  profile  data.  Using  this  method  we  lound  the  apparent  attenuation  for  vertical  component  L 
transverse  component  Lg ,  and  vertical  component  Pg  are  0(/>296/051,  C(/)^2()2/°”,  and  (?(/)=198/ 
respectively,  for  the  frequency  band  0.3</<10.0  Hz. 
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During  August  and  September,  1988,  the  United  Slates  and  the  Soviet  Union  conducted  a  Joint  Verification 
Experiment  in  which  each  countrv  was  permitted  to  observe  a  nuclear  explosion  at  the  other  country’s  test  site 
using  hydrodynamic  means.  We  have  scismically  recorded  the  Soviet  Explosion  at  four  stations  at  regional  dis¬ 
tance  (A<750  km).  The  pcak-to-pcak  particle  motion  in  the  P-  wave  varied  from  7.54x1 0-1  cm/scc  at  =  160  km 
to  1.02xl0-2  cm/scc  at  A=  740  km.  The  average  ntb (i.g )  was  5.97  corresponding  to  a  yield  of  120  KT,  based  on 
the  NTS  mb(Lt)  -  yield  relation  of  Patton[1988).  The  average  P„  comer  frequency  of  the  Soviet  explosion  was 
4  Hz  with  evidence  of  a  second  comer  at  about  0.7  Hz.  both  of  these  comers  could  be  the  effect  of  spall. 
Near-regional  seismograms  of  the  Soviet  JVE  explosion  show  a  large  arrival  within  the  Lt  phase  on  the 
transverse  component  at  all  four  sites.  This  may  correspond  to  tectonic  release  accompanying  the  explosion,  or 
to  the  effects  of  anisotropy  on  seismic  wave  propagation. 


Observations  of  10-,  20-,  and  30-1  lz  P-wave  spectral  amplitudes  from  earthquakes  and  explosions  arc  compared 
with  the  Archambcau  [1968,  1972]  earthquake  model  featuring  a  P-wave  falloff  of  to'3  beyond  the  comer  fre¬ 
quency,  a  modified  Brune  [1970,  1971]  earthquake  model  with  to'2  falloff,  and  the  Sharpe  [1942]  explosion 
model  which  has  a  to-2  falloff.  The  Archambcau  and  Shaqic  models  have  been,  in  part,  the  basis  of  a  proposal 
by  Evemden  et  al  [  19861  that  high-frequency  (=  30  Hz)  seismic  energy  could  provide  an  effective  solution  to  the 
problem  of  detection  and  identification  of  low-yield  coupled  and  fully  decoupled  underground  nuclear  explo¬ 
sions.  The  observations  of  earthquakes  show  an  increase  in  spectral  amplitude  with  moment  approximately  in 
agreement  with  the  to'2  falloff  model,  and  for  larger  moments,  in  disagreement  with  die  to'3  model.  Comparison 
of  theoretical  and  actual  seismograms  narrow-band  filtered  at  30-Hz  shows  that  in  part  the  increase  in  spectral 
amplitude  of  earthquakes  is  due  to  the  complex  and  long  duration  of  the  rupture  process  and  not  because  of  an 
increase  in  an  impulsive  first  arrival  like  that  characteristic  of  an  explosion.  The  30-11/  amplitudes  for  explo¬ 
sions  show  much  scatter,  and  many  events  have  a  spectral  falloff  greater  than  die  predicted  by  the  Sharpe 
model.  Whether  this  is  due  entirely  to  attenuation  or  is  the  actual  source  spectrum  is  not  determined.  High 
stress  drop  earthquakes  arc  predicted  to  have  larger  spectral  amplitudes  than  the  Shaqic  model.  Thus  any 
discrimination  technique  using  high-frequency  P-wave  spectra  should  probably  take  into  account  differences  in 
pulse  shape  and  amplitude  in  the  time  domain. 

Broadband  receiver  functions  determined  from  lelcsci.smic  P  waveforms  al  two  seismic  stations  in  eastern 
Kazakh,  U.S.S.R.,  were  inverted  for  the  veruca!  velocity  structure  beneath  the  stations.  1'hc  detailed  broadband 
receiver  functions  arc  obtained  by  stacking  source -equalized  radial  components  of  tclcscismic  P-waveforms.  A 
time-domain  inversion  of  the  radial  receiver  function  is  used  to  determine  the  structure  assuming  a  crustal  model 
parameterized  by  flat-lying,  homogeneous  layers.  The  general  features  of  the  inversion  results  are:  a  complex 
shallow  crust,  velocities  less  than  6  km/s  in  the  upper  crust,  a  high  velocity  (6.9-7. 5  km/s)  lower  crust,  and  a 
Moho  that  varies  between  47  and  57  km  depth.  These  results  compare  favorably  with  a  composite  velocity 
model  from  Deep  Seismic  Sounding  data,  but  show  lower  velocities  in  the.  upper  crust,  and  higher  velocities  in 
the  lower  crust.  The  results  indicate  that  the  crust  in  this  region  is  relatively  uniform,  however  the  nature  and 
depth  of  the  Moho  changes  significantly  across  the  region  The  receiver  function  inversion  structure  for  the 
Kazakh  sites  is  similar  to  a  published  receiver  function  structure  for  a  site  on  the  Canadian  shield. 
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Attenuation  of  Regional  Seismic  Phases  in  the  Great  Basin 


David  E.  Chavez1  and  Keith  F.  Priestley 
Seismological  Laboratory,  Mackay  School  of  Mines 
University  of  Nevada  -  Reno,  Reno,  Nevada  89557 


ABSTRACT 

Measurement  of  the  spat  :.,;  ;  xay  of  spectral  amplitudes  of  the  higher-mode 
seismic  surface  wave  tuiin  /  indicates  that. in  the  Great  Basin,  the  apparent 
seismic  quality  factor,  Q,  is  a  function  of  frequency.  Analysis  of  recordings 
from  a  NTS  nuclear  explosion  made  along  a  300  km  profile  through  northwest 
Nevada  yield  the  function  Q{f)  =  206/'  66  over  the  band  0.3  <  /  <  10.0 
Hz.  Similar  analysis  using  numerous  recordings  of  earthquakes  along  paths 
which  provide  a  good  average  sampling  of  the  Great  Basin  crust  gives  Q{f)  = 
214(±15)/°  54:1:0  09,  for  0.3  <  /  <  5.0  Hz.  If  the  crustal  sampling  by  Lg  en¬ 
ergy  from  nuclear  explosions  (surface  sources)  is  primarily  in  the  shallow  crust, 
as  has  been  suggested,  then  these  results  indicate  a  greater  frequency  depen¬ 
dence  of  apparent  Q  there  than  at  depth.  To  better  estimate  the  errors  in  Q 
at  each  frequency,  we  have  used  a  moving  window  spectral  analysis  technique 
described  by  Butler  et  al  (1987)  to  analysis  the  nuclear  explosion  profile  data. 
Using  this  method  we  found  the  apparent  attenuation  for  vertical  component 
Lg,  transverse  component  Lg,  and  vertical  component  Pg  are  Q(f)  =  296/°  51 , 
Q(f)  =  202/°  57,  and  Q(f)  =  198/°  45  respectively,  for  the  frequency  band 
0.3  </<  10.0  Hz. 

INTRODUCTION 

In  the  past  15  years,  considerable  efforts  have  been  devoted  to  the  under¬ 
standing  of  regional  seismic  phases.  The  majority  of  this  work  has  concentrated 
on  understanding  the  propagation  characteristics  of  the  crustal  phase  Lg.  Much 
less  effort  has  been  expended  on  the  other  prominent  regional  seismic  phases  P„, 
Sr, ,  Pg,  and  Sg.  Knowledge  of  the  att  nuation  charact  -istics  of  regional  seismic 
phases  is  important  for  both  earthquake  source  mechanism  studies  and  for  the 
discrimination  and  yield  estimation  of  small  nuclear  explosions.  Lg  consists  of  a 
superposition  of  higher  mode  surface  waves  and  because  of  its  averaging  prop¬ 
erties  in  excitation  and  transmission,  Nuttli  (1986)  has  suggested  that  it  can  be 
used  to  accurately  estimate  the  yield  of  nuclear  explosions.  However,  as  Nuttli 
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has  demonstrated,  it  is  necessary  to  know  the  Lg  attenuation  to  within  309? 
along  the  propagation  path  in  order  to  estimate  the  yield  to  within  10%.  Most 
Ls  attenuation  studies  to  date  have  concentrated  on  the  frequency  band  below  2 
to  3  Hz.  High  frequency  seismic  waves  have  taken  on  added  importance  recently 
because  they  may  be  more  diagnostic  els  to  different  earthquake  source  models, 
and  because  of  suggestions  that  they  may  improve  the  capability  to  distinguish 
earthquakes  and  small  or  decoupled  nuclear  explosions  (Evernden  et  al,  1986). 
Only  recently  has  data  become  available  for  investigating  the  attenuation  of 
higher  frequency  regional  seismic  phases.  In  this  study  we  examine  the  attenu¬ 
ation  of  the  Lg  propagating  within  the  Great  Basin  of  western  North  America. 
We  then  make  preliminary  steps  towards  understanding  the  attenuation  of  the 
Ps  phase  propagating  in  the  same  region. 

Cheng  and  Mitchell  (1981)  examined  the  spectra  of  long-period  Ls  for  three 
paths  in  the  Basin  and  Range  province  and  estimated  a  constant  shear  wave  Q 
of  85.  Singh  ind  Herrmann  (1983)  measured  the  decay  of  coda  waves  in  this 
same  region,  and  modeled  frequency  dependent  Q  as 

Q(f)  =  QoP  (l) 

where  Q0  is  the  Q  at  1  Hz.  They  found  Q(f)  =  250(±50)/°  45,±0  os\  and  were 
able  to  relate  their  coda  Q  values  to  the  Lg  spatial  attenuation  coefficients, 
concluding  that  the  same  mechanism  was  responsible  for  both  coda  decay  and 
Lg  attenuation.  Singh  and  Herrmann's  result  agrees  with  that  of  Cheng  and 
Mitchell  if  we  take  the  latter’s  value  cf  85  to  be  the  Lg  Q  at  5  seconds  period,  the 
predominate  period  in  their  data.  Peseckis  and  Pomeroy  (1984)  examined  the 
Lg  coda  of  several  nevada  Test  Site  (NTS)  explosions  recorded  at  regional  short 
period  stations  and  found  an  average  relationship  of  Q(f)  =  264/°  3.  However, 
Nuttli  (1986)  examined  a  similar  data  set  and  obtained  Q{f)  =  139/°  6. 

All  of  the  above  studies  used  band  limited,  analogue  data  recorded  by  only 
a  few  stations  and  along  a  limited  number  of  paths.  Using  broad  band  digital 
recordings  of  earthquakes  along  a  variety  of  paths,  as  well  as  data  from  a  profile 
recording  of  a  nuclear  explosion,  we  have  attempted  to  refine  the  Lg  attenuation 
model  for  the  Great  Basin,  and  measure  Pg  attenuation. 


Data 

The  data  analyzed  were  acquire  from  two  sources.  The  first  data  set  consist 
of  digital  seismograms  of  1 1  regional  earthquakes  in  and  around  the  Great  Basin 
recorded  at  four  broad-band  seismograph  stations  operated  by  the  Lawrence 
Livermore  National  Laboratory  (LLNL).  Figure  1  shows  the  locations  of  the 
earthquakes,  seismic  stations,  and  the  travel  paths  between  them.  We  selected 
earthquakes  which  were  recorded  by  at  least  three  of  the  four  stations  and 
which  covered  a  large  range  in  epicentral  distance  (i.e..  no  two  stations  were  at 
a  similar  distance).  As  indicated  by  the  figure,  the  travel  paths  used  provide  a 
good  sampling  of  the  Great  Basin  crust .  The  second  data  set  examined  consist  of 
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digital  recordings  of  a  nuclear  explosion  detonated  at  NTS  on  12/5/85.  Prior  to 
that  event,  fourteen  temporary  sites  were  instrumented  along  a  profile  covering 
the  distance  range  200  to  500  km,  extending  northwest  from  NTS  (Fig.  1).  A 
record  section  of  the  vertical  component  seismograms  from  the  nuclear  explosion 
is  shown  in  Figure  2.  The  Lg  windows  analyzed  are  marked  on  each  seismogram. 
The  beginning  of  each  window  corresponds  to  a  group  speed  of  3.7  km/sec,  while 
the  end  corresponds  to  the  arrival  of  the  fundamental  mode  surface  wave  train 
Rg,  typically  at  around  2.7  km/sec.  The  Rg  arrival  times  were  determined  by 
examining  a  low-pass  filtered  version  of  the  record  section.  The  same  criteria 
was  used  in  determining  the  Lg  windows  in  the  earthquake  data. 


Analysis 

Chavez  and  Priestley(1986)  determined  frequency  dependent  attenuation  for 
the  vertical  component  Lg  phase  from  these  data.  They  modeled  the  observed 
log-log  Lg  spectra  as 

-4(/,  R)  =  £0  5  €xp(  Q(/) )  (2) 

where  A(f ,  R)  is  the  spectral  amplitude  observed  at  a  distance  R,  f  is  frequency, 
t  is  the  travel  time,  S(f)  is  the  source  term,  and  Q  is  the  quality  factor.  Since 
this  model  does  not  consider  scattering  or  radiation  pattern  effects  it  provides  a 
measure  of  the  apparent  rather  than  the  intrinsic  Q.  Radiation  pattern  effects 
should  be  minimal,  however,  since  Lg  consists  of  a  large  number  of  rays  sampling 
a  major  fraction  of  the  focal  sphere.  Lg  has  been  successfully  modelled  as  surface 
wave  (Knopoff  et  al,  1973;  Panza  and  Calcagnille,  1975),  so  we  have  assumed 
that  the  frequency  domain  geometrical  spreading  scales  with  the  square  root  of 
distance.  Taking  the  logarithm  (base  10)  of  both  sides  of  (2)  yields 

logi0A(f,  R)  +  logl0R°  5  =  lo9l0S(f)-  |  |  /  (3) 

which  is  the  equation  for  a  straight  line,  with  the  source  term  as  the  intercept 
and  the  Q  term  controlling  the  slope.  Fixing  /,  we  know  A,  R ,  and  t  for  each 
of  the  stations  and  we  solve  for  S  and  Q  using  least-squares.  By  looping  over 
all  frequencies  we  obtain  the  source  and  Q  spectra. 

Figure  3  is  a  plot  of  the  Q{f)  functions  obtained.  Light  lines  indicate  results 
from  the  earthquake  data,  the  heavy  line  is  for  the  explosion  data.  Only  those 
values  whose  associated  linear  correlation  coefficient  exceeds  0.5  are  included. 
A  least-squares  fit  to  the  explosion  data  in  Figure  3  is 

<?(/)  =  206/°  C60  3  <  f  <  IQ.OH :  (4) 

while  a  weighted  fit  of  the  earthquake  data  gives 

Q(f)  =  214(±15)/° SM±0  o9)0.3  <  /  <  5.0//.-  (5) 
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The  latter  function  was  obtained  by  fitting  the  average  of  the  earthquake  derived 
Q  functions,  weighted  by  the  standard  deviation  at  each  frequency.  We  might 
expect  the  error  bars  for  the  explosion  data  to  be  no  greater  than  those  given 
in  (5)  for  the  earthquake  data.  If  we  assume  that  our  signal  to  noise  criteria 
was  successful  in  removing  random  error  from  the  spectral  amplitudes,  then  the 
error  bars  may  be  considered  as  a  measure  of  the  lateral  heterogeneity  of  the 
attenuation  structure.  The  relatively  small  magnitude  of  these  errors  suggests 
that  the  apparent  Q  function  varies  little  within  the  Great  Basin. 

The  Q  function  obtained  using  the  explosion  data  is  slightly  different  from  . 
that  based  on  the  earthquake  data.  While  the  values  of  Q0  found  for  both  data 
sets  are  essentially  identical,  the  explosion  data  indicate  a  greater  frequency 
dependence.  As  a  consequence,  the  Q  from  the  earthquake  data  are  consistently 
lower  than  that  from  the  explosion  data  for  frequencies  above  ]  Hz. 

One  drawback  in  using  the  above  method  to  measure  Q  is  that  it  is  difficult 
to  estimate  the  error  in  Q  at  each  frequency.  In  order  to  do  so,  we  have  repro¬ 
cessed  our  data  using  a  moving  window  spectral  analysis  technique  described 
by  Butler  et  al  (1987).  The  data  are  windowed  for  many  different  group  speeds 
using  overlapping  Gaussian  functions  whose  widths  are  scaled  by  the  epicentral 
distance  in  such  a  manner  that  the  windows  are  identical  with  respect  to  group 
speed  at  each  distance.  The  amplitude  spectra  for  each  window  are  then  them¬ 
selves  Gaussian  windowed  at  intervals  of  0.25  Hz,  and  the  windowed  spectrum 
is  integrated  to  yield  amplitude  as  a  function  of  group  speed  and  frequency. 
Figure  4  is  an  example  of  results  of  this  procedure.  Representing  seismic  data 
in  this  manner  provides  a  quantitative  description  of  the  frequency  content  of 
the  signal  as  a  function  of  time.  The  phases  Ps  and  Lg  are  clearly  visible  as 
ridges  in  frequency-group  speed  space.  The  loss  of  high  frequency  energy  in  the 
coda  is  also  apparent. 

Once  these  data  are  corrected  for  geometrical  spreading  we  can  solve  for  Q 
as  described  above,  only  this  time  we  obtain  apparent  attenuation  as  a  function 
of  group  speed  as  well  as  frequency.  The  attenuation  for  a  particular  phase  is 
determined  by  taking  the  average  over  the  group  speed  window  for  that  phase. 
The  standard  deviations  obtained  at  each  frequency  are  then  a  measure  of  the 
error  in  attenuation  for  that  phase. 

Before  doing  this  type  of  analysis  it  is  necessary  to  first  address  the  question 
of  the  nature  of  geometrical  spreading  when  using  this  technique.  Ls  is  com¬ 
posed  of  higher  mode  surface  waves  which  decay  in  the  frequency  domain  as 
R°  5.  However,  the  moving  window  spectral  analysis  method  employs  narrow 
time  windows  which  individually  do  not  contain  the  entire  phase.  If  any  energy 
at  a  particular  frequency  remains  outside  a  given  window,  then  the  amplitude 
at  that  frequency  and  group  speed  will  have  an  apparent  geometrical  spreading 
term  greater  than  0.5. 

We  have  investigated  the  effect  of  moving  window  spectral  analysis  on  ge¬ 
ometrical  spreading  with  an  experiment  using  synthetic  seismograms.  Lg  was 
synthesized  at  distances  and  azimuths  corresponding  to  the  stations  along  the 
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profile  by  summing  the  fundamental  and  first  25  Rayleigh  wave  modes  for  an 
attenuation-free  Great  Basin  model  (Priestley  and  Brune,  1978).  The  seismo¬ 
grams  were  then  processed  in  the  manner  described  above,  and  the  amplitude 
versus  distance  data  were  used  to  determine  apparent  geometrical  spreading 
as  function  of  frequency  and  group  speed.  The  results  are  plotted  in  Figure  5 
which  shows  that  geometrical  spreading  is  fairly  constant  for  those  frequencies 
and  group  speeds  for  which  there  is  signal.  Figure  5  also  shows  the  average 
apparent  geometrical  spreading  at  each  frequency  in  the  group  speed  window 
3.7  to  3  0  km/sec.,  corresponding  to  Lg.  It  can  be  seen  that  the  decay  term  is, 
for  the  most  part,  greater  than  0.5,  although  there  is  some  scatter. 

The  above  experiment  indicates  that  a  geometrical  spreading  term  of  0.6  is 
more  appropriate  than  0.5  at  most  frequencies,  although  at  other  frequencies 
terms  in  the  range  0.4  to  0.6  are  obtained.  In  order  to  determine  how  significant 
the  scatter  is,  we  processed  the  data  recorded  along  the  profile  assuming  geo¬ 
metrical  spreading  terms  of  0.0  to  1.2  in  steps  of  0.1  and  solved  for  Q.  Average 
Q  at  each  frequency  was  obtained  for  Lg  windows  of  3.7  to  3.0  km/sec  and  Pg 
windows  of  6.0  to  4.5  km/sec.  The  logQ  versus  log /  data  were  fit  by  least 
squares  to  give,  for  each  spreading  term,  the  function  Q(f)  —  Q}n  ■  Figure  6 
summarizes  the  results.  The  frequency  dependence  terms  obtained  for  vertical 
and  transverse  Lg  vary  little  with  assumed  geometrical  spreading,  and  the  Q0 
term  becomes  stable  for  spreading  terms  above  0.6.  This  suggests  that  (a)  the 
apparent  Q  of  both  vertical  and  transverse  Lg  is  frequency  dependent,  and  (b) 
the  scatter  in  spreading  term  about  0.6  evident  in  Figure  5  should  not  introduce 
any  serious  error. 

Figure  6  also  shows  that  Pg  behaves  differently  than  with  respect  to 
assumed  geometrical  spreading.  The  frequency  dependence  term  for  this  phase 
consistently  decreases  with  increasing  spreading  while  the  Qa  term  increases. 
However,  if  we  assume  that  the  actual  Q  of  Pg  behaves  as  Qafn .  then  the  linear 
correlation  coefficients  shown  in  Figure  6c  indicate  that  an  apparent  geometrical 
spreading  of  less  than  0.7  is  required. 

Assuming  an  apparent  geometrical  spreading  term  of  0.6,  we  obtained  Q 
along  the  profile.  Figure  7  shows  the  results  for  each  frequency  and  group  speed 
and  averaged  over  windows  corresponding  to  Lg  and  Pg.  The  fits  to  the  data 
over  the  frequency  band  0.3  <  /  <  10.0  Hz  are 


Q(f)  =  296/°  51 

(6) 

Q{})  ~  202/° 57 

(7) 

Q(f)  =  198  f0Ai 

(8) 

for  vertical  component  Lg,  transverse  component  Lg,  and  vertical  component 
Pg,  respectively.  Vertical  Lg  has  Q  similar  to  both  vertical  Pg  and  transverse 
Lg  for  frequencies  below  about  1.5  Hz.  Above  that,  vertical  Lg  Q  is  consistently 
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greater  than  the  others,  except  above  10  Hz,  where  both  vertical  and  transverse 
Lg  again  have  similar  Q.  There  is  a  pronounced  bump  in  all  the  plots  at  2  to  3 
Hz.  Since  this  feature  appears  for  all  the  phases  considered,  we  suspect  that  it 
might  be  an  artifact  due  to  recording  site  conditions.  Finally,  vertical  Lg  and 
vertical  Pg  both  have  Q  independent  of  frequency  for  the  range  4  to  9  Hz. 


Discussion 

Our  results  are  in  reasonable  agreement  with  those  of  Cheng  and  Mitchell 
(1981)  and  Singh  and  Herrmann  (1983).  Both  (4)  and  (5)  yield  values  of  Q 
which  are  within  25%  of  that  found  by  Cheng  and  Mitchell,  if  we  assume  their 
results  to  reflect  the  Q  at  5  seconds.  Our  Q0  values  are  within  the  range  given  by 
Singh  and  Herrmann  (1983).  The  frequency  dependence  obtained  here  using  the 
earthquake  data  is  in  very  good  agreement  with  theirs,  if  we  take  into  account 
the  uncertainties  in  the  two  values.  Again,  the  explosion  data  yield  a  Q  function 
w-ith  a  greater  frequency  dependence. 

Our  findings  do  not  support  the  conclusions  of  either  Peseckis  and  Pomeroy 
(1984)  or  Nuttli  (1986).  Peseckis  and  Pomeroy's  Qc  is  greater  than  ours,  and 
we  find  no  evidence  to  suggest,  a  frequency  dependence  as  low  as  f°  3.  Nuttli's 
result  of  f°  6  is  closer  to  ours:  however,  his  Q0  is  too  low,  according  to  our 
data.  Perhaps  these  discrepancies  are  indicative  of  the  larger  scatter  inherent 
in  time  domain  measurements  of  frequency  dependent  quantities.  Since  all  of 
the  previous  studies  considered  only  a  few  paths,  it  is  likely  that  the  results 
presented  here  are  more  representative  of  the  true  apparent  Q  in  the  Great 
Basin. 

The  discrepancy  between  the  frequency  dependence  inferred  from  the  explo¬ 
sion  versus  the  earthquake  data,  if  significant,  suggest  that  energy  from  the  two 
source  types  does  not  sample  the  crust  in  an  identical  manner.  Campillo  et  al 
(1985)  present  ray  tracing  diagrams  which  show  that  the  crustal  sampling  by 
Lg  waves  is  influenced  by  source  depth.  Their  results  suggest  that  Lg  energy 
from  surface  sources  (like  explosions)  propagates  primarily  in  the  shallow  layers, 
while  an  earthquake  at  10  km  depth  generates  energy  which  uniformly  samples 
the  entire  crust.  If  this  is  the  case  in  the  great  Basin,  then  our  data  indicate 
that  there  is  greater  frequency  dependence  of  Q  in  the  shallow  crust  than  at 
depth. 
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Fig.  1.  Map  showing  the  region  sampled  by  our  data.  The  permanent  Lawrence 
Livermore  stations  are  indicated  by  large  squares,  dots  show  the  epicenters  of 
the  earthquakes  recorded  at  those  stations  and  considered  in  this  study.  Small 
squares  are  the  temporary  sites  instrumented  to  record  the  December  5,  1985 
NTS  nuclear  explosion. 


Fig.  2.  Reduced  record  section  for  the  vertical  component  data  recorded  along 
the  profile  in  Figure  1.  The  windows  used  to  compute  the  Ls  spectra  are  marked 
on  each  record.  The  change  in  signal  character  at  490  km  is  due  to  the  different 
instrument  response  there. 


Fig.  3.  Great  Basin  Q(f)  functions.  The  light  lines  are  the  results  using 
earthquake  data,  the  heavy  line  is  from  the  explosion  profile. 

Fig.  4.  Examples  of  moving  window  spectra  for  the  closest  (top)  and  furthest 
(bottom)  stations  along  the  nuclear  explosion  profile  in  Figure  1.  Contours  are 
log  ground  amplitude  in  microns. 


Fig.  5.  Apparent  geometrical  spreading  along  the  profile  determined  by  invert¬ 
ing  synthetic  data.  Bottom  figure  is  average  apparent  geometrical  spreading  for 
the  indicated  group  speed  window.  The  dashed  line  is  at  0.5,  the  solid  line  at 
0.61  is  the  average. 


Fig.  6.  Effect  of  geometrical  spreading  on  determination  of  Q(f)  by  applying  the 
moving  window  spectra  technique  to  the  data  recorded  along  the  profile.  Closed 
circles  denote  vertical  component  Lg,  squares  denote  transverse  component  Lg, 
and  triangles  denote  vertical  component  Pg.  Arrows  in  (a)  and  (b)  indicate 
the  term  solved  for  by  fitting  average  log  Q  and  log  frequency  values  for  group 
speed  windows  appropriate  for  the  phases  Lg  and  Pg\  (c)  is  the  linear  correlation 
coefficient  obtained  from  the  least-squares  fits. 


Fig.  7.  Apparent  Q  functions  for  vertical  and  transverse  Lg  and  vertical  Ps. 
obtained  by  the  moving  window  spectral  analysis  technique  assuming  a  geo¬ 
metrical  spreading  term  of  0.6.  Error  bars  represent  the  standard  deviation  in 
apparent  Q  at  each  frequency  for  averages  over  3.7  to  3.0  km/sec  for  Lg  and 
6  0  to  4.5  km/sec  for  Pg.  Least-squares  fit  to  the  data  yield  the  relationships 
indicated. 
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Abstract.  During  August  and  September,  1988,  the  United 
States  and  the  Soviet  Union  conducted  a  Joint  Verification 
Experiment  in  which  each  country  was  permitted  to  observe 
a  nuclear  explosion  at  the  other  country’s  test  site  using 
hydrodynamic  means.  We  have  seismically  recorded  the 
Soviet  explosion  at  four  stations  at  regional  distances 
(A  <,  750km),  and  here  report  those  observations.  The 
peak-to-peak  particle  motion  in  the  P-wave  varied  from 
7.54x10*'  cm/sec  at  A  =  160  km  to  1.02xl(TJ  cm/s ec  at 
A  =  740  km.  The  average  mb(Ls)  was  5.97  corresponding 
to  a  yield  of  1 1 8  KT,  based  on  the  NTS  m^g)  -  yield  rela¬ 
tion  of  Patton  [1988],  The  average  Pn  comer  frequency  of 
the  Soviet  explosion  was  4  Hz  with  evidence  of  a  second 
comer  at  about  0.7  Hz.  Both  of  these  comers  could  be 
affected  by  spall.  Near-regional  seismograms  of  the  Soviet 
JVE  explosion  show  a  large  arrival  within  the  L,  phase  on 
the  transverse  component  at  all  four  sites.  This  may 
correspond  to  tectonic  release  accompanying  the  explosion, 
or  to  the  effects  of  anisotropy  on  seismic  wave  propagation. 

Introduction 

During  August  and  September,  1988,  the  United  States 
and  the  Soviet  Union  conducted  a  Joint  Verification  Experi¬ 
ment  (JVE)  for  the  purpose  of  calibrating  yield  measure¬ 
ments  at  their  nuclear  weapons  testing  sites.  Scientists  from 
each  country  were  permitted  to  observe  a  nuclear  explosion 
at  the  other  country’s  test  site  using  hydrodynamic  means. 
In  addition,  as  pan  of  the  official  experiment  each  country 
made  teleseismic  (A  >  2500  km)  measurements  at  five 
national  seismograph  stations  in  the  U.S.  and  and  five  in  the 
Soviet  Union.  The  Natural  Resources  Defense  Council 
-  NRDC)  reached  an  agreement  with  the  Soviet  Academy  of 
■’■>'.  it  cite  tS/s.'>)  to  make  addinonal  seismic  recordings  of  the 
:>->ict  A'£  explosion  in  the  Soviet  Union  at  regional 
ci-stance  ranges  (<  2500  km).  Here  we  present  preliminary 
analysis  of  the  seismic  data  recorded  at  four  sites  in  Eastern 
Kazakhstan  at  distances  less  than  750  km. 

Seismic  Instrumentation 

Seismic  data  for  the  Soviet  JVE  nuclear  explosion  were 
recorded  at  four  sites  in  Eastern  Kazakhstan.  Three  seismo- 
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graphs  were  at  the  previously  occupied  NRDC-SAS  sites 
[Berger  et  al,  1987]  surrounding  the  Soviet  test  site  at  near- 
regional  distance  range:  Karasu  (KSU,  A=  160  km),  Kai 
karlinsk  (KKL,  A  =  255  km),  and  Bayanual  (BAY,  A  =  255 
km).  The  fourth  instrument  was  located  at  the  SAS  seisr.x 
graph  site  at  Talgar  (TLG,  A  =  740  km)  to  the  south  of  the 
nuclear  weapons  testing  site.  The  locations  of  the  recording 
sites  and  of  the  explosion  are  shown  in  Figure  1. 

The  four  seismograph  sites  lie  within  the  Kazakh  fold  sys¬ 
tem,  a  belt  of  folded  and  faulted  Paleozoic  rocks  [Leith, 
1987],  The  region  was  marked  by  extensive  igneous 
activity  during  the  late  Hercynian  Orogeny  in  the  Permo- 
Triassic  Period.  KSU,  KKL,  and  BAY  are  located  on  simi 
lar  granitic  intrusions  of  Paleozoic  to  early  Mesozoic  rocks. 
The  granites  were  intruded  into  deformed  Paleozoic  sedi 
mentary  rocks.  Subsequent  tectonic  activity  has  resulted  in 
regional  crustal  uplift  and  thickening.  The  SAS  seismo¬ 
graph  station  at  TLG  is  situated  on  the  northern  flanks  of 
the  Zaili-Alatau  mountains  and  is  sited  within  Precambnan 
and  Lower  Paleozoic  crystalline  rocks.  The  most  prominent 
structural  feature  in  the  area  is  the  Chingiz  fault  which  cuts 
northwest  to  southeast  across  the  Soviet  test  site.  There  are 
indications  of  rejuvenated  movement  along  the  fault  during 
the  Neogene  and  Quaternary  [Zhuravlev  and  Uspensky, 
1971],  Geological  exposures  in  the  vicinity  of  the  Soviet 


Fig.  1.  Location  of  the  four  regional  seismograph  stations 
KKL,  BAY,  KSU,  and  TLG,  and  the  location  of  t’-r 
Soviet  JVE  explosion. 
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test  site  suggest  the  area  is  fairly  heterogeneous.  Despite 
the  geological  complexity,  Ruzashin  et  al.  [1977]  find 
efficient  Lg  propagation  across  the  Kazakh  fold  system. 

The  characteristics  of  the  three  near-regional  sites  KSU, 
KKL,  and  BAY  have  been  described  by  Berger  et  al. 
[1988].  The  permanent  NRDC-SAS  seismographs  previously 
operating  at  KKL,  BAY,  and  KSU  [Berger  et  al,  1987; 
1988]  had  been  removed  in  late  1987.  We  reinstalled  short- 
and  long-period  seismographs  in  each  of  the  vaults  to  record 
the  Soviet  JVE  nuclear  explosion.  The  short-period  seismo¬ 
graphs  consisted  of  three  component  Teledyne  Geotech  S-13 
1  second  free  period  seismometers  and  EDA  PRS-4  digital 
data  loggers.  The  long-period  seismographs  consisted  of  a 
three  component  Teledyne  Geotech  SL2 10/220  15  second 
free  period  seismometers  and  a  Kinemetrics  PDR-2  digital 
data  logger.  All  these  data  were  sampled  at  200  Hz,  with 
the  exception  of  the  short-period  seismograph  at  BAY 
which  was  sampled  at  100  Hz.  At  TLG  the  sensors  are 
SM-3,  1.6  sec  free  period  seismometers.  Data  are  recorded 
at  100  Hz  on  the  PUSK-2,  a  120  dB  uynamic  range  digital 
data  logger  developed  jointly  by  the  Soviet  and  Georgian 
Academy  of  Sciences. 


Data 


Examples  of  the  short-period  vertical  component  seismo¬ 
grams  for  the  Soviet  JVE  explosion  are  shown  in  Figure  2. 
The  KKL  seismogram  shows  an  emergent,  low  frequency 
first  arrival,  followed  about  three-quarters  of  a  second  later 


2  57X10-3  ~-2  13X1  °" 

— - \l 

KKL 


4- 


2  56X10-J 
— vW' 

BAY 


C  l U  - .  . 

- mm " !  . . i\v-'iV — ■ — 


— V-Www\fvVw^ 


KSU 


1  59X10-’  ■*  1  02X10 

— V  o'w<ll  I 

TGL  _ 


Fig.  2.  Shon-period  vertical  seismograms  of  the  Soviet 
JVE  explosion  recorded  at  the  four  regional  seismograph 
sites.  The  upper  part  of  each  section  of  the  figure  shows 
10  seconds  of  the  P-wave,  the  lower  portion  shows  100 
seconds  of  the  whole  wave  train  at  KKL,  BAY  and  KSU 
and  150  seconds  at  TLG.  Two  expanded  plots  of  the 
KSU  P-wave  are  shown,  the  upper  is  the  original  data, 
and  the  lower  is  the  same  data  section  after  low-pass 
filtering  with  a  butterworth  4  pole  filter  with  the  comer  at 
1.5  Hz.  The  peak-to-peak  particle  motion  for  initial 
arrival  and  for  the  largest  amplitude  in  the  first  ten 
seconds  is  noted  and  has  units  of  criVsec. 


by  a  low  amplitude,  higher  frequency  arrival.  KKL  is  past 
the  Pn  -  Pg  cross-over  distance  predicted  from  the  eastern 
Kazakhstan  velocity  structure  [Antonova  et  al.,  1978;  and 
Priestley  et  al.,  1988],  and  we  assume  the  low  frequency 
first  arrival  is  Pn.  The  Pg-wave  train  is  dominated  by  the 
large  amplitude  dilatational  arrival  on  the  vertical  and  radial 
components  following  the  first  Pn  arrival  by  about  1.5  sec. 
The  dilatational  motion  of  this  arrival  suggest  a  supercritical 
reflection  from  the  Moho.  Synthetic  reflectivity  seismo¬ 
grams  computed  for  the  Priestley  et  al.  [1988]  velocity 
structure  agree  in  both  the  Pn-PmP  interval  time  ana  in  the 
relative  amplitudes  observed  at  KKL,  as  well  as  for  obser¬ 
vations  of  the  same  phase  at  BAY  and  KSU.  There  is  a 
well  developed  Lg  phase  seen  on  the  vertical  component, 
but  as  discussed  below,  the  Lg  amplitudes  are  largest  on  the 
transverse  component. 

The  BAY  seismogram  shows  a  more  impulsive,  larger 
amplitude  and  higher  frequency  Pn  arrival  than  that 
observed  at  KKL.  However,  like  KKL,  the  BAY  P-wave  is 
dominated  by  a  large  amplitude  dilitationa)  arrival  about  1.5 
seconds  after  the  initial  P-wave  arrival  (see  Figure  4).  The 
vertical  component  Lg  at  BAY  is  similar  to  that  seen  a: 
KKL.  The  maximum  amplitude  in  the  Lg  phase  observed  at 
BAY  is  also  on  the  transverse  component.  However,  most 
of  the  Lg  energy  in  the  transverse  BAY  record  occurs  in  a 
single  pulse. 

The  closest  recording  of  the  Soviet  JVE  was  made  at 
KSU,  within  the  Pn-Pg  crossover  distance.  Both  the  short 
and  long-penod  recordings  from  KSU  show  the  characteris¬ 
tic  3  Hz  resonance  associated  with  the  site  (upper  KSU 
recording  in  Figure  2).  Neither  previous  borehole 
recordings  [Berger  et  al.,  1988],  nor  nearby  recordings 
(=  200m)  outside  the  vault  [Myers,  personal  communication, 
1988]  show  as  strong  a  resonance  indicating  that  it  is  not 
characteristic  of  the  site  in  general,  but  of  the  pier  in  the 
vault.  We  have  also  included  a  lowpass  filtered  time  series 
to  remove  the  effect  of  the  site  resonance.  This  filtered 
seismogram  of  KSU  shows  a  large  amplitude  arrival  follow¬ 
ing  the  first  arrival  by  about  2  seconds,  corresponding  to  the 
time  for  the  PmP  reflection  at  this  distance. 

The  TLG  recording  was  made  at  a  much  larger  distance 
than  the  other  recordings.  The  largest  phase  on  the  vertical 
seismogam  has  an  amplitude  of  4.27  x  10'2  cm/sec,  with  a 
group  velocity  of  7.18  km/s.  There  is  a  well  developed  Lg 
phase  at  TLG  with  the  maximum  amplitude  in  the  Lg  wave 
train  on  the  transverse  component. 

Preliminary  analysis 

Figure  3  shows  the  instrument  corrected  (solid  line), 
whole-path  Q-corrected  (dashed  line),  and  noise  (dotted 
line)  spectra  for  the  Pn,  Pg,  and  Lg  phases  from  the  long- 
period  vertical  component  seismogram  recorded  at  BAY. 
The  KKL  spectra  are  almost  identical  to  those  computed 
from  the  BAY  recording.  The  dashed  Lg  spectra  is 
corrected  for  Q  using  our  unpublished  attenuation  results 
(Q  l,  =  367  f  °48).  We  have  made  the  often  used  assump¬ 
tion  (e.g.  Evemden  et  al.  [1986])  that  Q  Pj  /  Q  l,  =  2.25,  so 
Q  p(  =  825  f  °'4*,  in  plotting  the  dashed  Pg  spectrum.  The 
dashed  Pn  spectrum  is  corrected  for  attenuation  using  the 
value  of  Q  =  9000  suggested  by  Evemden  et  al.  [1986). 
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Fig.  3.  Pn,  P{,  and  Lg  spectra  for  the  vertical  component 
seismogram  recorded  at  BAY.  The  time  windows  used 
are  indicated  above  the  time  series.  A  10  percent  cosine 
taper  was  applied  prior  to  transforming.  The  comer  fre¬ 
quency  is  noted  on  the  spectra.  The  solid  curve  is  the 
spectrum  of  the  indicated  phase.  The  dashed  curve  is  the 
attenuation  corrected  spectrum  using  the  values  quoted  in 
the  text  The  dotted  curve  is  the  spectrum  of  the  noise. 
For  Pn  and  P,  the  noise  spectra  was  computed  by 
transforming  an  equal  length  time  series  immediately 
before  the  first  P-wave  arrival;  for  Lg  an  equal  length 
time  series  immediately  before  the  L,  arrival  was  used. 


The  Pn  and  Pg  spectra  shows  a  good  signal  to  noise  ratio  to 
at  least  50  Hz.  If  we  parameterize  the  Pn  and  Pg  spectra 
with  a  two  line  model  between  1  and  30  Hz,  the 
unconecied  spectra  show  similar  behavior  with  comer  fre¬ 
quencies  of  about  4  Hz  and  falloff  slopes  above  the  comer 
approximately  proportional  to  to"3.  There  is  a  suggestion  of 
a  second  comer  in  the  spectra  at  about  0.7  Hz.  The 
corrected  spectra  show  about  the  same  comer  frequency  and 
a  slightly  shallower  high  frequency  falloff.  If  the  actual 
attenuation  for  Pg  is  closer  to  Q  as  some  scattering  stu¬ 
dies  sugg:st,  the  falloff  slope  becomes  closer  to  of2.  The 
Lg  spectrum  shows  signal  above  the  noise  to  about  3  Hz. 

The  average  comer  frequency  for  the  Soviet  explosion  is 
somewhat  higher  than  comer  frequencies  of  similar  size 
nuclear  explosions  detonated  at  NTS  recorded  at  comparable 
distances.  The  lower  comer  frequencies  of  the  U.S. 
explosions  may  reflect  the  the  higher  crustal  and  upper  man¬ 
tle  attenuation  in  the  Basin  and  Range,  or  possibly  a  lower 
coupling  coefficient  of  the  weaker  rocks  in  the  Nevada  Test 
Site.  However,  for  NTS  explosions  of  about  the  same  size, 
spall  affects  frequencies  in  the  range  0.5  to  2.0  Hz  [Taylor 
and  Randall,  1988].  We  have  no  knowledge  of  the  spall 
characteristics  of  the  Soviet  test  site  and  cannot  eliminate 
differences  in  spall  between  the  two  testing  areas  as  the 
cause  of  the  differences  in  the  spectra. 


Figure  4  shows  the  transverse  component  seismograms 
after  rotating  the  horizontal  components  to  radial  and 
transverse.  These  seismograms  show  significant  energy  at 
the  Sg  /  Lg  arrival  time.  The  BAY  seismogram  shows  a 
remarkably  strong  coherent  SH  arrival  on  the  transverse 
component.  This  is  the  largest  amplitude  arrival  on  the 
BAY  seismogram.  KKL  and  KSU  also  show  significant  SH 
energy,  however  this  is  spread  out  over  several  seconds 
duration.  The  filtered  BAY  seismograms  show  a  distinct 
Love  wave  arrival  on  the  transverse  component  preceding 
the  Rayleigh  wave  on  the  vertical  and  radial  components. 
Previous  studies  have  noted  anomalous  surface  wave  radia¬ 
tion  from  Semipaliunsk  nuclear  explosions  [Rygg,  1979] 
These  observations  have  been  attributed  to  tectonic  release 
equivalent  to  motion  on  thrust  faults  dipping  45°.  However, 
there  have  been  suggestions  that  apparent  quadrupole 
seismic  radiation  from  explosions  may  arise  from  the  aniso 
tropic  properties  of  the  source  or  propagation  medium  (e  g. 
Mandal  and  Toksoz,  1989). 

Nuttli  [1973]  defined  a  magnitude  mb(Lg)  based  on  the 
amplitudes  of  1  second  Lg  waves,  and  Nuttli  [1986a]  argued 
for  the  usefulness  of  mb(Lg)  in  estimating  yields  of  NTS 
nuclear  explosions.  Nuttli  [1986b,  1987]  used  Lg  arrivals 
from  Soviet  nuclear  explosions  in  eastern  Kazakhstan  which 
were  recorded  in  western  Europe  and  southern  Asia  to  esti¬ 
mate  yields  of  Soviet  explosions.  The  greatest  uncertainty 
in  the  mb  -  yield  relationship  arises  from  uncertainties  in  the 
attenuation  function  used  for  correcting  the  Lg  amplitudes 
for  the  effects  of  anelasticity.  Since  we  have  made  direct 
measurements  of  the  Lg  attenuation  only  in  the  region  of  the 
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Fig.  4.  Transverse  component  recordings  of  the  Soviet 
JVE  at  KSU,  KKL  and  BAY  are  shown  at  the  lop  of  the 
figure.  The  KSU  trace  is  the  short-period  record  low  pass 
filtered  at  1.5  Hz.  The  remaining  traces  are  long-period 
broadband  seismograms.  Also  shown  for  comparison  are 
the  radial  and  vertical  components  at  BAY.  At  the  bot¬ 
tom  of  the  figure  low-pass  filtered  (comer  of  0.2  Hz) 
transverse,  radial  and  vertical  BAY  seismograms  aie 
shown.  Each  trace  is  self-scaled.  The  filtered  BAY 
seismograms  show  significant  Love  waves  The  absolute 
amplitude  of  the  filtered  BAY  transverse  component  is 
about  four  times  smaller  than  the  vertical  component. 
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Soviet  test  site,  we  have  used  only  Lg  recordings  from 
BAY,  KKL  and  KSU  in  estimating  mj,(L4).  The  for 

the  Soviet  JVE  is  5.68910.002.  Nuttli  (1986a,  1986b,  1987) 
argues  for  the  portability  of  -  yield  relationships  derived 
from  NTS  data;  however,  Patton  [1988]  cautions  that  cou¬ 
pling  differences  may  have  a  significant  influence.  An 
mb(Lj)  of  5.689  corresponds  to  a  yield  of  148  KT  using  the 
Nuttli  [1986]  mb(Lf)  -  yield  relation  and  118  KT  using  the 
Patton  [1988]  mb(Lg)  vs.  yield  relation,  both  of  which  are 
based  on  NTS  data.  Sykes  and  Ekstrom  [1989]  using  a  com¬ 
bined  mb  and  M,  magnitude  -  yield  relation  have  estimated 
the  yield  of  the  Soviet  JVE  to  be  113  KT. 

Summary 

Broadband  regional  seismograms  of  the  1988  Soviet  JVE  at 
four  sites  in  the  U.S.S.R.  have  been  shown.  The  P  wave 
spectra  show  good  signal-above-noise  up  to  at  least  50  Hz. 
The  transverse  components  show  significant  energy  at  the 
Sg  /  Lg  arrival  time,  and  Love  waves  at  longer  periods.  This 
Soviet  JVE  data  provides  an  opportunity  to  test  methods  of 
analysis,  developed  in  other  areas,  on  data  recorded  wholly 
within  the  U.S.S.R.  Using  one  such  method,  the  mb(Lg)  - 
yield  relations  developed  for  NTS,  we  have  estimated  the 
yield  of  this  event. 

Acknowledgements.  We  thank  Steve  Taylor,  George 
Zandt,  and  George  Randell  for  helpful  discussions.  Christina 
Hackman  of  the  University  of  Colorado  and  Kenneth  Smith 
of  the  University  of  Nevada  helped  in  the  data  collection 
effort.  The  data  collection  was  supported  by  a  gram  from 
the  Natural  Resources  Defense  Council.  This  research  was 
partially  supported  by  grants  from  the  Air  Force  Geophysi¬ 
cal  Laboratory  under  contract  number  F19628-87-K-009,  the 
Lawrence  Livermore  National  Laboratory  under  contract 
number  W-7405-ENG-48.  and  the  National  Science  Founda¬ 
tion  under  contract  number  EAR8708506. 

References 

Antonova,  L.  V.,  F.  F.  Aptikayev,  R.  1.  Kurochkina,  I.  L. 
Nersesov,  A.  V.  Nikolayev,  A.  I.  Ruzaykin,  Y.  N. 
Sedova,  A.  V.  Sitnikov,  F.  S.Tregub,  L.  D.Fedorskaya, 
and  V.  1.  Khalturin,  Experimental  Seismic  Investigation  of 
the  Earth's  Interior,  AS  USSR,  Institute  of  Physics  of  the 
Earth,  Publishing  House  "Nauka",  Moscow,  155p,  1978. 
Berger,  J.,  J.  N.  Brune,  P.  A.  Bodin,  J.  S.  Gomberg,  D.  M. 
Carrel,  K.  F.  Priestley,  D.  E.  Chavez,  W.  R.  Waller,  C.  B. 
Archambeau,  T.  B.  Cochran,  L.  L.  Nersesov,  M.  B. 
Gokhberg,  O.  A.  Stolyrov,  S.  K.  Daragen,  N.  D.  Taras- 
sov,  Y.  A.  Sutelov,  A  new  US  -  USSR  seismologica!  pro¬ 
gram.  EOS,  68,  105,  130-111,  3987. 

Berger,  J.,  H.  K.  Eissler,  F.  L.  Vernon,  I.  L.  Nersesov,  M. 
B.  Gokhberg,  O.  A.  Stolyrov,  and  N.  T.  Tarasov,  Studies 
of  high-frequency  seismic  noise  in  Eastern  Kazakhstan, 
Bull.  Seismo.,  Soc.  Am.,  78,  1744-1758,  1988. 


Evemden,  J.,  C.  Archambeau,  and  E.  Cranswick,  An  evalua¬ 
tion  of  seismic  decoupling  and  underground  nuclear  test 
monitoring  using  high  frequency  seismic  data,  Rev  Geo- 
phys.,  24,  143-215,  1986. 

Leith,  W.,  Geology  of  NR DC  seismic  stations  in  Eastern 
Kazakhstan,  USSR,  USGS  Open-File  Report  87-597, 
1987. 

Mandal,  B.,  and  M.  Toksoz,  Radiation  patterns  from  explo¬ 
sions  in  anisotropic  media,  paper  presented  at  DOE/LLNL 
symposium  on  explosion- source  phenomenology,  March 
14-16,  1989,  Lake  Tahoe  California. 

Nuttli,  O.,  Seismic  wave  attenuation  and  magnitude  relations 
for  eastern  North  America,  J.  Geophys  Res.,  78,  876-888, 
1973. 

Nuttli,  O.,  Yield  estimates  of  Nevada  Test  Site  explosions 
obtained  from  seismic  Lg  waves,  J.  Geophys.  Res  ,  91, 
2737-2151,  1986a. 

Nuttli,  O.,  Lg  magnitudes  of  selected  East  Kazakhstan 
underground  explosions.  Bull.  Seism  Soc.  Am.,  76,  1241- 
1251,  1986b. 

Nuttli,  O.,  Lg  magnitudes  at  Dcglen,  East  Kazakhstan  under¬ 
ground  explosions,  Bull.  Seism.  Soc.  Am.,  77,  679-681, 
1987. 

Patton,  H.  J.,  Application  of  Nuttli’s  method  to  estimate 
yield  of  Nevada  Test  Site  explosions  recorded  on 
Lawrence  Livermore  National  Laboratory's  digital  seismic 
system,  Bull.  Seism.  Soc.  Am.,  78,  1759-1772,  1988. 

Priestley,  K.  F.,  G.  Zandt,  and  G.  E.  Randall,  Crustal  struc¬ 
ture  in  Eastern  Kazakh,  U.S.S.R.  from  teleseismic 
receiver  functions,  Geophys.  Res.  Lett.,  15,  613-616, 
3988. 

Rygg,  E„  Anomalous  surface  waves  from  underground 
explosions,  Bull.  Seism.  Soc.  Am.,  69,  1995-2002,  1979. 

Ruzkin,  A.  I.,  I.  L.  Nersesov,  V.  I.  Khaturin,  and  P.  Molnar, 
Propagation  of  Lg  and  lateral  variations  in  crustal  struc¬ 
ture  in  Asia,  J.  Geophys.  Res.,  82,  307-376,  1977. 

Sykes,  L.,  and  G.  Ekstrom.  Comparison  of  seismic  and 
hydrodynamic  yield  determinations  for  the  Soviet  joint 
verification  experiment  of  1988,  Proc.  Natl.  Acad  Sci. 
USA,  86,  3456-3460,  1989. 

Taylor,  S.  R.,  and  G.  E.  Randall,  The  effects  of  spall  on 
regional  seismograms,  Geopys.  Res.  Lett.,  16,  211-214, 
1989. 

Zhuralev,  B.  Y.,  and  E.  P.  Uspensky,  Large  faults  of  the 
Chingiz  range,  American  Geophysical  Union  translation. 

Geotectonics,  5,  305-309,  1971. 


K.  Priestley  and  W.  Walter,  Seismologica)  Laboratory, 
University  of  Nevada,  Reno,  Reno,  NV  89557. 

V.  Martynov  and  M.  Rozhkov,  Institute  of  the  Physics  of 
the  Earth,  Soviet  Academy  of  Science,  Moscow,  U.S.S.R. 


(Received:  August  25,  1989; 
Accepted:  September  27,  1989) 


20 


•  HtSICAL  RtSLARCI I.  VOL  93,  NO  86.  1‘ACF.S  6318-Lj.N,  JUNE  10,  ',981' 


Jk  >l  ; 


Observations  of  High-Frequency  P  Wave  Earthquake  and  Explosion  Spectra 
Compared  With  to-3,  to-2,  and  Sharpe  Source  Models 

Wiluam  R.  Walter,  James  N.  Brune,  and  Keith  F.  Priestley 

Seismological  Laboratory,  Mackay  School  of  Mines,  University  of  Nevada  at  Reno 


Jon  Fletcher 

U  S.  Geological  Survey,  Menlo  Park,  California 


Observations  of  10-,  20-,  and  30-Hz  P  wave  spectral  amplitudes  from  canhquakcs  and  explosions  are  com¬ 
pared  with  the  Archambeau  [1968,  1972)  earthquake  model  featuring  a  P  wave  falloff  of  to'3  beyond  the  comer 
frequency,  a  modified  Brune  [1970,  1971]  earthquake  model  with  to'2  falloff,  anti  the  Sharpe  [1942]  explosion 
model  which  has  a  to-2  falloff.  The  Archambeau  and  Sharpe  models  have  been,  in  part,  the  basis  of  a  proposal 
by  Evernden  el  al.  [1986]  that  high-frequency  (=  30  Hz)  seismic  energy  could  provide  an  effective  solution  to 
the  problem  of  detection  and  identification  of  low-yield  coupled  and  fully  decoupled  underground  nuclear  ex¬ 
plosions.  The  observations  of  earthquakes  show  an  increase  in  spectral  amplitude  with  moment  approximately 
in  agreement  with  the  to-2  falloff  model  and,  for  larger  moments,  in  disagreement  with  the  to-3  model.  Com¬ 
parison  of  theoretical  and  actual  seismograms  narrow-band  filtered  at  30  Hz,  shows  that  in  part  the  increase  in 
spectral  amplitude  of  earthquakes  is  due  to  the  complex  and  long  duration  of  the  rupture  process  and  not  be¬ 
cause  of  an  increase  in  an  impulsive  first  arrival  like  that  characteristic  of  an  explosion.  The  30-llz  amplitudes 
for  explosions  show  much  scatter,  and  many  events  have  a  spectral  falloff  greater  than  the  or2  predicted  by  the 
Sharpe  model.  Whether  this  is  due  entirely  to  attenuation  or  is  the  actual  source  spectrum  is  not  determined. 
High  stress  drop  earthquakes  are  predicted  to  have  larger  spectral  amplitudes  than  the  Sharpe  model.  Thus  any 
discrimination  technique  using  high-frequency  P  wave  spectra  should  probably  take  into  account  differences  in 
pulse  shape  and  amplitude  in  the  lime  domain. 


Introduction 

Thir  paper  compares  observations  of  P  wave  spectral  ampli¬ 
tudes  at  10,  20,  and  30  Hz  from  numerous  earthquakes  and  explo¬ 
sions,  covering  a  range  of  magnitudes,  with  several  theoretical 
predictions.  In  particular  the  spectral  amplitudes  arc  compared 
with  the  Sharpe  [1942]  explosion  model  and  with  the  Archam¬ 
beau  [1968,  1972]  earthquake  stress  relaxation  model.  The  differ¬ 
ence  in  spectral  amplitudes  between  these  two  models  is  the  bads 
of  a  method  for  discriminating  between  earthquakes  and  under¬ 
ground  nuclear  explosions  suggested  by  Evernden  el  al.  [1986],  A 
preliminary  report  on  the  results  of  these  studies  was  given  by 
Brune  et  al.  [1986a], 

Evernden  et  al.  [1986]  proposed  that  an  effective  solution  to  the 
problem  of  detection  and  identification  of  low-yield  coupled  and 
fully  decoupled  underground  explosions  was  available  using 
l’:gh  frequency  seismic  data  ranging  up  to  30  or  40  Hz.  The  pro¬ 
posal  is  based  on  the  theoretical  predictions  of  radiated  P  wave  se¬ 
ismic  spectra  corresponding  to  the  Archambeau  [1968,  1972] 
model  for  earthquake  sources  and  the  Sharpe  [  1942]  model  for  ex¬ 
plosion  sources.  The  Archambeau  [1968,  1972]  theoretical  model 
fot  an  expanding  stress-relaxation  (failure)  zone  has  been  refined 
and  developed  by  Minster  [1973],  Stevens  [1980],  and  Archam¬ 
beau  and  Minster  [1978]  among  others.  Various  models  used 
nonsphcrically  expanding  or  translating  spherical  failure  zones  to 
simulate  time-dependent  ellipsoidal  failure  zones.  Stevens  [1980] 
showed  that  the  volume  relaxation  source  formulation  used  in  ear¬ 
lier  work  could  be  expressed  in  terms  of  failure  boundary  surface 
integrals  (stress  pulse  equivalents)  Thus  in  the  limit  it  seems  that 
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these  volume  stress-relaxation  models  could  be  modified  to  give 
the  same  results  as  the  more  physically  reasonable  models  involv¬ 
ing  stress  relaxation  on  a  plane,  c.g.,  as  in  finite  clement  modeling 
of  Archuleta  and  Frazier  [1978],  or  idealized  approximate  stress 
pulse  models  such  as  that  of  Brune  [1970,  1971],  However,  when 
one  takes  into  account  the  inherent  complexities  in  the  actual  run 
turn  process,  it  is  unlikely  that  any  of  these  idealized  models  ade¬ 
quately  represents  the  spectra  of  most  earthquakes.  Much  evi¬ 
dence  indicates  that  earthquake  ruptures  are  complex  multiple 
events  [IFyss  and  Brune,  1967;  Trifunac  and  Brune,  1970; 
Kartamori  and  Stewart,  1978].  Recognition  of  this  fact  has  stimu¬ 
lated  numerous  theoretical  and  numerical  models  involving  com¬ 
plexity  [Haskell,  1964;  Aki,  1972;  Brune,  1970;  Blandford,  1975; 
Papageorgiou  and  Aki,  1983a,  b\  Kanamori  1979;  Munguia  and 
Brune,  1984].  Recent-high  quality  data  from  the  Anza  seismic  ar¬ 
ray  [ Berger  et  al.,  1984;  Fletcher  et  al.,  1987]  have  indicated  r- r>* 
this  complexity  extends  down  to  small  events  [Brune  el  al.  19b;  l 
Frankel  et  al.,  1986]. 

The  theoretical  P  wave  explosion  spectra  used  by  Evernden  et 
al.  [1986]  is  that  proposed  by  Sharpe  [1942J  and  independently 
developed  by  Latter  [1961].  In  the  Sharpe  theory  the  explosion  is 
modeled  as  an  arbitrary  pressure  function  applied  to  the  interior  of 
a  spherical  surface  at  an  "clastic  radius"  from  the  source.  The 
elastic  radius  is  defined  as  the  distance  from  the  detonation  point 
to  the  surface  beyond  which  the  medium  responds  linearly  (elasti 
cally)  and  thus  is  a  purely  medium-dependent  property  for  a  given 
maximum  pressure  amplitude.  For  a  perfectly  decoupled  or  ove- 
decoupled  explosion  in  a  cavity,  the  clastic  radius  is  the  cavity  i 
dius.  For  a  tamped  or  underdccouplcd  explosion,  the  elastic  rarihs 
is  larger  than  the  perfectly  decoupled  cavity  radius  for  the  same 
yield  explosion,  since  the  coupling  factor  in  a  solid  medium  is 
larger  than  lire  coupling  factor  in  air.  The  Sharpe  model  as  t  •  i 
by  Evernden  et  al.  ]  1986]  (and  as  developed  by  Latter  ]  1961 ))  fins 
a  step  or  Heaviside  function  in  pressure  acting  at  the  clastic  ra- 
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dins.  Olhcr  clastic  radius  models  such  as  that  of  Mueller  and  Mur¬ 
phy  [1971]  use  a  pressure  function  with  an  initial  spike  which  de¬ 
cays  exponentially  to  a  constant  value.  This  model  has  more 
high-frequency  energy  al  and  above  the  comer  frequency  than  the 
pressure  step  function  Sharpe  model.  Evernden  el  al.  [1986]  argue 
based  on  the  Mueiler-Murphy  model  and  finite  difference  calcula¬ 
tions  by  Patterson  [1966]  that  the  pressure  step  function  Sharpe 
model  is  a  conservative  estimate  of  the  high-frequency  P  wave  en¬ 
ergy  radiated  by  both  lamped  and  decoupled  explosions.  We  will 
coniine  our  analysis  in  this  paper  to  tamped  explosions. 

One  uf  the  factors  of  crucial  importance  in  the  Evcmdcn  et  al. 
proposal  is  the  shape  of  the  P  wave  spectrum  beyond  the  corner 
frequency.  The  Archarnhcau  [1968,  1972 1  model  for  earthquakes 
(with  the  rupture  velocity  approximately  equal  to  the  shear  wave 
velocity)  predicts  a  high-frequency  falloff,  beyond  the  comer  fre¬ 
quency  pro|x>rlional  to  to'3.  For  simplicity  in  this  paper  we  will 
refer  to  die  model  as  the  W3P  mode!  (to-3  falloff  for  P  waves). 
This  kind  of  falloff  has  llie  paratloxic.il  result  dial  for  a  given 
stress  drop,  very  large  cardiquakcs  do  not  radiate  any  more  high 
frequency  energy  than  small  earthquakes  ( mb~  1.8).  This  is  a 
consequence  of  die  fact  that  die  low  frequency  spectral  level  is 
proportional  to  the  cube  of  the  source  dimension,  which  is  exactly 
cancelled  out  by  the  oT3  falloff  from  the  comer  frequency  (which 
scales  inversely  with  the  source  dimension).  Evernden  el  al. 

1 1986]  propose  that  the  co“3  falloff  of  the  P  wave  spectrum  for  the 
Archambeau  (1968,  1972J  model  for  cardiquakcs,  being  different 
from  the  co  2  falloff  if  the  Sharpe  model,  enhances  the  discrimina¬ 
tion  achievable  from  the  differences  in  P  wave  comer  frequencies 
between  compatible  cardiquakcs  and  explosions  (compatible  in 
that  they  have  die  same  low-frequency  P  wave  amplitudes). 

It  should  be  noted  that  die  W3P  model  as  used  by  Evernden  cl 
al.  does  not  lake  into  account  phase  complexity  of  earthquakes 
per  se,  since  only  the  spectrum  is  used  and  not  the  actual  time 
domain  signal  (or  phase).  Thus  a  conclusion  regarding  the  ampli¬ 
tude  spectra  of  the  W3P  model  for  earthquakes  docs  not  necessari¬ 
ly  argue  one  way  or  die  olhcr  as  to  whether  time  domain  measure¬ 
ments  might  be  used  to  distinguish  between  earthquakes  and  ex¬ 
plosions,  since  earthquakes  are  complex  long-duration  ruptures, 
whereas  explosions  arc  impulsive  stress  pulses  on  the  inside  of 
some  actual  or  equivalent  cavity.  (The  interrelations  between 
spectra,  peak  amplitude,  complexity,  and  duration  of  the  signal, 
and  the  high  frequency  falloff  have  been  discussed  by  Hanks 
[  1979|.) 

The  problem  of  determining  the  slope  of  the  spectral  falloff  at 
high  frequencies  is  greatly  complicated  by  attenuation,  scattering, 
and  sile  effects.  The  determination  of  seismic  moment,  being 
based  on  much  lower  frequencies,  is  relatively  insensitive  to  these 
effects.  Thus  if  data  from  one  source  region  with  a  constant  set  of 
stations  and  ray  palhs  arc  compared,  any  increase  in  high- 
frequency  spectral  amplitude  with  moment  can  be  attributed  to  the 
source.  Therefore  this  paper  will  address  the  variation  of  10-,  20-. 
and  30  Hr.  earthquake  and  explosion  P  wave  spectral  amplitudes 
with  seismic  moment  al  10- km  hyixiceritral  distances  without 
making  any  correction  for  Q. 


Tiii:ori:tk:ai.Curvi-;s 

In  the  Evernden  cl  al.  (1986|  paper.  Figure  14  shows  the 
tlieorelie.il  P  wave  spectra  of  the  W3P  model.  The  theoretical 
spectra  have  a  constant  amplitude  up  to  the  corner  frequency  ff, 
given  hy  [ Minster ,  197.7] 

/■f  =  — !—( 3(i2vr)''3  (1) 

2  7t  L 


and  then  falls  off  as  to  3.  Here  L  is  the  maximum  rupture  dimen¬ 
sion,  vr  is  die  rupture  velocity,  and  a  is  the  P  wave  velocity.  The 
relation  between  the  moment  and  the  low-frequency  asymptote  for 
the  W3P  model  was  derived  from  the  low-frequency  (point 
source)  approximation  for  a  double  couple  dislocation  [Kalis- 
Dor  ok,  1959J 
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and  the  moment  expression  for  the  W3P  model  (assuming  a  Pois¬ 
son  solid)  given  by  | Minster ,  1973;  Minster  and Suteau,  1977] 

..  60 n  , ,  ... 
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In  (2)  p  is  the  medium  density,  R  is  the  hypocentral  distance,  RM 
is  the  radiation  pattern,  M 0  is  the  seismic  moment,  |)  is  the  S  wave 
velocity,  and  III  and  Uj  are  die  S  and  P  wave  zero  frequency 
spectral  amplitudes.  In  (.3)  o  is  the  stress  drop  (or  die  preslrcss 
level).  Equaling  (2)  and  (3)  and  solving  for  ijJ  gives 
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where  p  is  the  rigidity.  This  result  was  checked  with  Figure  22 
from  Evernden  cl  al.  1 1986),  which  shows  earthquake  spectra 
predicted  by  die  W3P  model  and  explosion  spectra  predicted  by 
the  Sharpe  model.  The  vertical  axis  in  their  Figure  22  gives  rela¬ 
tive  amplitude  on  the  left-hand  scale  and  mb  on  the  right-hand 
scale.  To  relate  mb  with  moment  and  low-frequency  s]x:clral 
amplitude  we  have  used  the  mb  versus  the  12,,  spectral  level  rela¬ 
tionship  of  K.  F.  Priestley  and  J.  N.  Brunc  (Spectral  scaling  for  the 
Mammoth  Lakes,  California,  earthquakes,  submitted  to  die  Geo¬ 
physical  Journal  of  the  Royal  Astronomical  Society ;  hereafter 
refered  to  as  PB88).  The  values  obtained  agreed  with  those 
predicted  by  (4)  within  the  uncertainty  of  determining  mb.  This 
was  further  checked  with  die  graph  of  mb  versus  A/0  given  by 
Tucker  and  Prune  [1977]  for  San  Fernando  aftershocks,  with 
almost  identical  results. 

The  values  used  for  the  W3P  model  arc  a  -  6.0  ktn/s,  p  =  2.8 
g/cm3,  p  =  3.36  x  10"  dyn/cm,  and  a  Poisson  medium  is 
assumed.  We  have  further  assumed  that  die  rupture  velocity  is 
equal  to  the  shear  wave  velocity.  As  (1)  indicates,  decreasing  the 
rupture  velocity  causes  the  comer  frequency  to  decrease. 
Evernden  el  al.  1 1986)  show  in  their  Figure  16  dial  as  the  rupture 
velocity  decreases  below  the  shear  wave  velocity,  the  spectrum 
develops  an  intermediate  falloff  between  die  constant  low- 
frequency  value  and  the  to  3  falloff  al  frequencies  much  higher 
dian  the  comer  frequency.  Thus  the  case  of  die  rupture  velocity 
equal  to  die  shear  wave  velocity  gives  die  maximum  s|vclral 
amplitudes  for  P  waves  at  all  Irequencics  predicted  by  the  W3P 
model  for  a  given  stress  drop  and  moment.  (For  rupture  rates 
greater  than  the  shear  wave  velocity,  the  Archambeau  ( 1968, 
1972)  model  predicts  a  falloff  proportional  to  ok2  beyond  the 
comer  frequency.) 

For  the  to"2  model  of  earthquake  spectra,  we  used  the  Brunc 
[1970,  1971]  model  of  source  spectra  with  die  added  assumption 
dial  the  P  wave  to  .9  wave  amplitude  ratio  al  low  frequencies  is 
(a/|3)3.  The  ratio  is  appropriate  for  die  low-frequency  (point 
source)  approximation  for  a  double  couple  dislocation  given  tn 
(2).  We  have  also  assumed  a  comer  frequency  shift  for  P  waves 
relative  to  S  waves  of  i,  as  indicated  in  several  studies  [e  g  . 
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.  . ;  i  lie  ■  c  ,'s  .!  foe  the  medium  arc  tlie  same  as  the 

W3P  model.  Hctealtcr  for  simplicity  this  model  will  be  referred 
to  .is  i!>  ’  W2I1  model  (to  2  falloff  for  P  waves). 

The  S'  i'pe  model  curves  represent  both  lamped  and  decoupled 
e\ j  f  ision.;  T  be  moment  used  for  the  Sharpe  curves  was  derived 
from  the  low  frequency-moment  relation  of  (2),  as  if  the  source 
were  an  earthquake,  using  the  same  parameters  for  the  recording 
medium  as  the  W3P  model.  In  this  paper  we  have  used  a  max¬ 
imum  pressure  of  elastic  response  of  130  bars  as  assumed  by 
Fi'cmu/ <  •!  i  t  it!  |  !<>K6|  in  estimating  the  spectral  amplitude  from 
explosions  at  the  Nevada  test  sue.  We  have  also  assumed  shot 
medium  parameters  of  (X  =  3.5  km/s  and  p  =  2.0  g/em3.  All  the 
theoretical  curves  presented  in  this  pajvr  have  been  nonnalized  to 
10-km  l:\poccnlial  distance  and  multiplied  by  a  factor  of  2  to 
approximate  the  free  surface  cflect. 

1).\1  a 

Theie  are  numerous  sources  of  good  high-frequency  data  now 
available  as  a  result  of  the  advent  of  high  dynamic  range  digital 
event  recorders  [e.g.,  Brune  cl  al.,  19X0)  and  telemetered  arrays 
1  Berger  it  nl „  19X4;  Fletcher  el  al.,  1987],  The  earthquake  data 
used  in  this  study  are  recordings  from  the  permanent  Anza  array 
in  southern  California  [Berger  et  al.,  1984;  Fletcher  et  al.,  1987], 
and  tempirary  digital  instruments  deployed  to  record  aftershocks 
of  the  iv7X  Oaxaca,  Mexico,  earthquake  [Munguia  cl  al.,  1978) 
and  aftershocks  of  the  1980  Mammoth  Lakes  earthquakes 
| Archuleta  et  al.,  1982;  Priestley  et  al.,  1985;  PB88).  In  order  to 
compare  the  data  with  the  source  models  with  as  little  distortion 
due  to  attenuation  and  site  effects  as  possible  we  have  used  data 
recorded  on  solid  rock  with  hypocenlral  distances  as  short  as  pos¬ 
sible.  All  the  data  have  been  normalized  to  a  10-kin  hypocentra! 
distance  and  have  not  been  corrected  for  Q. 

The  data  presented  here  are  all  far-field  data  in  the  frequency 
domain,  except  for  the  1985  Nliehoacan  main  shock.  The  instru¬ 
ments  used  in  this  study  have  a  flat  response  in  velocity  or 
acceleration  between  5  and  35  Hz.  The  spectra  were  dctemiincd 
from  vertical  component  recordings.  The  data  lime  scries  range 
from  3/4  s  to  several  seconds  after  die  initial  P  arrival  and  ending 
before  the  .S'  wave  (except  for  the  Michoucan  event).  The  time 
series  were  tapered  with  a  tO'v  cosine  taper  at  each  end.  A  fast 
F  ourier  transform  was  applied  to  the  tapered  time  scries;  the  mag¬ 
nitude  of  the  result  was  dien  rotated  (dividing  hy  at  for  velocity 
and  cr  for  acceleration)  to  displacement.  The  Anza  spectra  were 
then  smoothed,  and  the  spectra!  amplitude  determined  directly. 

’  i  r  '!i  r  ■  .  of  the  data,  lit-'  magnitude  of  the  spectral  amplitude  at 

i'-s’b  i.'.."  f  vq.iency  was  delennmed  by  averaging  about  a  6-Hz 
window  i  entered  on  the  frequency  ol  interest  (he.,  27-33  Hz  for 
T>  H ■)  using  unsmoothed  spectra.  Where  more  than  one  station 
recording  was  available  for  each  event,  we  followed  Archuleta  el 
al  |  14X21  by  taking  the  average  of  the  log  of  the  spectral  ampli¬ 
tude,  rather  than  a  simple  arithmetic  average. 

I  he  At;/. i  array  is  a  II)  statu  a,  three -component,  16-bil  digi- 
tallv  telemetered  array,  deployed  along  the  San  Jacinto  fault  zone 
in  southern  California.  We  have  analyzed  events  with  seismic 
moments  ranging  from  lO17  to  about  1()21  dyn  cm  and  w  ith  hypo- 
central  distances  from  5  to  20  km.  The  Anza  array  sites  arc  on 
'olid  rock,  and  we  cxjiect  little  distortion  due  to  attenuation  al 
these  dict.inces.  Hough  et  al.  (S.  E.  Hough,  J.  G.  Anderson,  J. 
Brune,  1  Vernon.  I.  Berger,  J.  I  lelcher.  L  Harr,  T.  Hanks  and  L. 
Maker.  An-  mi.ition  near  An/a.  Cahloinia.  submitted  to  Bulletin  of 
t’u  -SVn/c,. ‘git  al  Si’i  .( ly  i‘f  Ami  rim)  and  S.  E.  Hough  and  J.  G. 
Aralei  cii  iljhscrc ations  of  high  (requeue y  spec u a  al  An/a,  Cali¬ 


fornia:  Implications  lor  Q  stniclure,  submitted  t  '  P.tS.i  .■  ,h. 
Scismological  Society  of  America)  have  found  high  (9  value  -  lor 
Anza  sites.  Recent  comparison  of  spectra  recorded  in  a  K-‘  1 

deep  borehole  and  at  the  surface  by  llaar  et  d.  |!9X7)  • 
amplifications  up  to  a  factor  of  5  in  the  10-  to  30  H/  • 
(beyond  the  expected  factor  of  2  at  a  free  surface)  and  significant 
attenuation  of  frequencies  above  40  H/.  However,  this  is 
expected  to  affect  all  events  in  a  similar  manner  and  should  not 
change  comparisons  of  events  recorded  al  the  same  stations.  We 
also  have  strong  motion  data  for  two  earthquakes,  an  M  ~  S  O  a, 
an  M  =  4.7,  that  occurred  prior  to  the  installation  of  the  Ae/.i 
array. 

The  Oaxaca  aftershock  data  were  recorded  on  an  eight-soil 'on 
digital  array  set  up  after  the  main  shock  of  November  29,  h'/.v 
[Reyes  cl  al.,  1978).  These  sites  were  also  located  on  solid  uuA  o 
minimize  attenuation  and  site  effects.  These  events  ate  g  - 
deeper  titan  events  at  An/.a  since  they  arc  associated  with  tae  sub 
duelion  zone.  Tltey  range  in  hypocenlral  distances  from  it)  m  if. 
km. 

The  Mammoth  events  were  recorded  on  digital  seismographs 
deployed  during  the  Mammoth  Lakes  earthquake  sequence  in 
1980  [ Archuleta  el  al.,  1982).  To  minimize  attenuation  we  used 
events  recorded  at  the  three  stations  CON,  FIS,  and  MGE  l.n  ,u,d 
on  rock.  The  hypocenlral  distances  of  these  events  range  from  5  . 
20  km. 

Wc  have  used  one  very  large  event,  'die  Michoacan,  Mexico 
eartliquake  of  September  19,  1985.  A  network  of  strong  motion 
accclcrographs  located  on  hard  rock  sites,  was  instill'  -.!  in  the 
region  prior  to  the  event  [Anderson  ct  al.,  1986).  Three  o!  the 
accclcrographs  were  directly  over  the  aftershock  zone  and  inferred 
rupture.  We  have  used  the  spectra  of  the  entire  vertical  component 
of  these  stations.  As  a  first-order  approximation  to  normalize  to 


TABLE  1.  Parameters  for  Pahutc  Mesa,  Nevada,  Test  Site  Explosions  used 
in  This  Study 


Event  Name 

Depth 

m 

Yield* 

kt 

mb 

Hypocenlral 
distance  km 

Station 

Alcmcndro 

1064 

566 

6  1 

3,2 

L  '34 

3.2 

L-C4 

3.4 

Echo  Pc-‘ 

5.1 

L-05 

6.1 

L-02 

10.0 

L-06 

12.6 

L-01 

Halfbcak 

819 

300 

6.1 

2.1 

S -5 

2.1 

S-7 

Ink-. 

siy 

290 

6.0 

1.6 

O  X  . 

3.3 

S 

Scotch 

977 

150 

5.7 

6.1 

6-0* 

Boxcar 

1166 

1200 

6.3 

10.4 

S-034 

22.5 

S-0,4 

Farm 

689 

75 

5.5 

4.2 

4 

4.7 

6 

64 

9.7 

T 

21  1 

W  o 

21.2 

W-.  1 

Col  wick 

633 

75 

5.5 

18  3 

W-| 

23.4 

W-!8 

207 

W  " 

Pcpato 

680 

75 

5.5 

22.5 

V'- 

IS  4 

A  <> 

Shccpshcad 

640 

100 

5.6 

15  1 

\V- 

17.7 

W  1  ’ 

22.8 

w-i  • 

*  Yields  from  Springer  and  Kinnman  119711.  Barker  el  al.  ( 1  uft  -  i  - 
estimated  from  mb. 


23 


Walter  et  al.:  Observations  of  High-Frequency  P  Waves 


6321 


10-km  hypocenlral  distance,  we  have  assumed  the  energy  radiated 
per  unit  area  is  constant  at  the  fault  surface,  and  that  energy  is 
conserved.  Therefore  we  multiply  the  spectral  amplitude  by  die 
square  root  of  the  ratio:  twice  die  fault  surface  area  divided  by  die 
area  of  a  sphere  at  10  km.  This  correction  increases  the  spectral 
amplitude  by  about  a  factor  of  3.7. 

The  explosion  data  used  here  arc  underground  nuclear  explo¬ 
sions  from  Pahutc  Mesa  on  the  Nevada  test  site  (NTS).  We  have 
several  explosions  ranging  in  size  from  approximately  100  to 
1 200  kilotons  (kl),  recorded  on  velocity  and  acceleration  sensors. 
The  event  names,  approximate  yields,  and  hypocentral  distances 
are  given  in  Table  1.  The  hypocentral  distances  range  from  3  to  25 
km.  The  moment  for  each  explosion  was  determined  using  the 
formula 

A70  =  18.135 +  mt  (5) 

taken  from  Priestley  and  Brunc  (PB88)  (assuming  the  Qq  -  M 0 
relation  of  (2)  and  die  parameters  of  du  W3P  model).  The  velo¬ 
city  recordings  at  Paliutc  Mesa  are  surface  recordings  on  sedi¬ 
ments.  Using  the  velocity  model  of  Leonard  and  Johnson  [1987) 
and  data  from  the  shot  medium  (D.  L.  Springer,  personal  com¬ 
munication  1987),  we  expect  an  impedance  amplification  of  about 
1 .4-1.6  for  events  at  650-m  depth  and  about  1.8  for  1 100-m  depth 
(based  on  die  difference  between  (pot)1'2  at  die  surface  and  at 
depth).  In  addition  the  attenuation  is  expected  to  be  greater  for 
these  events  than  for  the  earthquake  data.  To  obtain  an  estimate 
for  the  attenuation  of  the  P  wave  explosion  data  we  assume  a  (fre¬ 
quency  independent)  value  of  Q  =  50  above  600-m  depdi  and  Q  = 
100  for  depths  greater  than  600  m  (L.  Johnson,  personal  commun¬ 
ication,  1987).  Events  recorded  at  5.0-km  cpiccntral  distances 
have  travel  times  of  about  1.6  s.  Assuming  a  travel  lime  of  about 
1.2  s  in  die  Q  =  100  material  and  0.4  s  in  the  Q  =  50  material 
gives  a  total  amplitude  decrease  of  about  6.6  at  30  Hz  and  1.9  at 
10  Hz.  There  is  evidence  diat  at  high  frequencies,  attenuation  is 
frequency  dependent;  this  could  be  expected  to  decrease  the 
attenuation  al  30  Hz  relative  to  that  at  10  Hz.  However,  the 
attenuation  in  this  region  is  not  vet  well  determined,  so  we  have 
not  attempted  to  correct  for  either  attenuation  or  sediment 
amplification.  Based  on  the  assumption  for  Q  cited  above,  we 
expect  the  uncorrecled  spectral  amplitudes  to  approximate  die 
source  sjiectra  al  10  Hz  and  underestimate  it  at  30  Hz. 

For  events  Colwiek,  Pepalo,  Shccpshcad,  and  Farm  we  have 
accelerograph  recordings.  These  acceleration  recordings  arc  either 
oil  hard  rock  surface  sites  (stations  W-i 7  and  W-18)  or  from  500- 
m-deep  boreholes  (station  W-9).  The  spectral  amplitude  from  the 
borehole  station  was  multiplied  by  a  factor  of  2  to  approximate  a 
surface  recording.  It  is  interesting  to  note  dial  at  die  more  distant 
(hypocentral  distances  from  25  to  33  km)  crystalline  rock  surface 
and  borehole  statioas  the  30  Hz  P  wave  energy  was  lost  in  the 
noise.  For  the  event  Farm,  we  have  bodi  hard  rock  accelerograph 
recordings  (hypocentral  distances  approximately  21  km)  and  four 
accelerograph  recordings  on  sediments  (hypocentral  distances  4-9 
kill).  The  S|icclral  amplitudes  of  the  recordings  on  sediments  are 
approximately  a  factor  of  2  larger  at  10,  20,  and  30  Hz  dian  the 
more  distant  recordings  on  rock.  These  recordings  on  rock  may 
underestimate  die  source  amplitude  due  to  impedance.  All  of  die 
explosion  data  have  been  normalized  to  10-km  hypocentral  dis¬ 
tance. 

Analysis 

The  difference  in  high-frequency  behavior  of  die  W2P,  W3P, 
and  Shar|x:  source  models  as  the  event  size  increases,  is  illustrated 
in  Figure  1,  where  spectral  amplitude  at  10,  20,  and  30  Hz  is  plot¬ 


ted  versus  seismic  moment.  The  Sharpe  1 1 942 )  model  is 
represented  by  the  dashed  curve  and  was  calculated  using  a  max¬ 
imum  pressure  of  clastic  resjxinsc  of  130  bars.  The  W2P  and  W3P 
models  are  represented  by  solid  lines.  The  stress  drop  in  bars  asso¬ 
ciated  with  each  curve  is  indicated  by  o.  The  W3P  curve  was 
derived  from  (1)  and  (4)  assuming  the  rupture  velocity  is  equal  to 
die  shear  wave  velocity.  The  data  have  all  been  normalized  to  10- 
km  hypocentral  distance  and,  as  discussed  in  the  introduction,  not 
corrected  for  Q. 

The  theoretical  curves  all  show  spectral  amplitude  increasing 
linearly  with  moment  up  to  about  1015-102U  dyn  cm.  This  is 
expected  for  any  source  model,  since  for  these  small-moment 
events  the  wavelength  at  die  frequency  of  interest  (10,  20,  or  30 
Hz)  is  large  compared  with  the  source  dimension.  In  other  words, 
diese  small  events  have  comer  frequencies  higher  dian  die  fre¬ 
quency  being  plotted,  so  the  spectral  amplitude  scales  directly 
with  moment.  For  the  larger-moment  (above  about  1020  dyn  cm) 
events  the  W3P  model  predicts  the  spectral  amplitudes  will 
remain  constant  for  a  constant  stress  drop.  The  W2P  model  in 
contrast  predicts  for  diese  large-moment  events  that  the  spectral 
amplitudes  will  increase  with  moment  approximately  as  A/J'5  for 
constant  stress  drop.  The  Sharpe  curve  has  similar  scaling  to  the 
W2P  model. 

Figure  1  a  shows  10-Hz  S|>ectral  amplitude  versus  moment.  The 
Oaxaca  events  show  more  10-11/  energy  than  An/a  events  of 
comparable  moment  and  dins  could  be  interpreted  as  having 
higher  stress  drops  than  the  An/a  events  in  either  the  W2P  or 
W3P  models.  However,  the  two  models  give  different  absolute 
values  for  stiess  drops.  For  example,  the  largest  moment  An/a 
events  would  be  interpreted  as  having  a  stress  drop  of  about  10 
bars  in  the  W2P  model,  while  the  W3P  model  would  give  stress 
drops  of  about  400  bars.  Many  of  die  Oaxaca  events  have  as  much 
or  more  10-Hz  energy  dian  the  Sharpe  curve.  The  explosion 
points  scatter  about  die  Sharpe  curve  with  die  recordings  on  rock 
having  less  10-Hz  energy  than  predicted  by  the  model.  The  six 
points  representing  die  five  small,  75-150  kt,  events  (Farm  was 
recorded  both  on  rock  and  sediments)  have  about  an  order  of 
magnitude  range  in  10-Hz  amplitude. 

As  we  look  al  higher  frequencies  (20  Hz  in  Figure  16  and  30 
Hz  in  Figure  lc)  the  difference  in  the  predicted  spectral  amplitude 
of  the  W2P  and  W3P  models  for  large-moment  events  becomes 
more  pronounced.  The  earthquake  data  (Anza,  Oaxaca,  Mam¬ 
moth,  and  Michoacan)  show  a  clear  increase  in  spectral  amplitude 
with  moment,  approximately  in  agreement  widi  the  W2P  model, 
and  in  disagreement  with  a  constant  stress  drop  W3P  model  for 
die  larger  events.  If  we  were  to  interpret  die  data  with  the  W3P 
model,  stress  drop  would  have  to  increase  significantly  with 
moment.  For  the  largest  moment  events  such  as  the  Michoacan 
main  shock  (Figure  lc)  the  W3P  model  implies  a  stress  drop  of 
the  order  of  a  kilobar.  This  is  in  distinct  disagreement  w  ith  static 
.stress,  apparent  stress,  and  dynamic  stress  as  calculated  by  Ander 
son  et  al.  [1986]  and  Priestley  arul  Masters  [1986]  of  about  5-20 
bars. 

The  NTS  explosions  show  increasing  scatter  as  we  look  at  die 
higher  frequencies  (Figures  16  and  lc).  The  hard  rock  recordings 
and  die  Farm  accelerograph  recording  on  sediments  have  lower 
spectral  amplitudes  at  30  Hz  (Figure  lc)  by  1-2  orders  of  magni¬ 
tude.  One  possible  explanation  for  this  is  source  depth.  The  higher 
spectral  amplitude  events  arc  deeper  (800-1 100  in)  than  the  lower 
spectral  amplitude  events  (630-690  m).  The  average  depdi  to  the 
water  table  in  Pahutc  Mesa  is  alxnil  650  m  \Lconard  ami  Johnson, 
1987).  Events  below  the  water  table  might  be  expected  to  couple 
more  high-frequenev  energy  into  seismic  waves.  The  majority  o! 
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Fig.  1  a 

Fig.  1.  V  wave  displacement  spectral  observations  of  earthquakes  (solid 
symbols)  and  NTS  explosions  plotted  as  a  function  of  moment,  and  com¬ 
pared  with  theoretical  curves.  All  symbols  and  curves  normalized  to  a  hy- 
poccntral  distance  of  10  km.  No  attenuation  correction  was  applied  to  the 
data,  Events  where  more  than  one  station  recording  was  used  have  rela¬ 
tively  larger  symbol  size.  Theoretical  curves  given  for  constant  stress  drop 
denoted  by  ft  in  bars.  W3P  model  given  for  constant  rupture  velocity  equal 
to  the  shear  wave  velocity.  The  Sharpe  curve  was  calculated  assuming  a 
maximum  pressure  of  130  bars  and  is  indicated  by  the  dashed  curve,  (a) 
The  10-Ilz  P  wave  displacement  versus  moment,  (6)  20-Hz  P  wave  dis¬ 
placement  spectral  amplitude  versus  moment,  and  (c)  30-Hz  P  wave  dis¬ 
placement  tpcctral  amplitude  versus  moment.  Note  the  increase  in  spec¬ 
tral  amplitude  with  moment  approximately  in  agreement  with  a  constant 
stress  drop  W2P  ( oT2  falloff  beyond  the  comer  frequency)  model  and  for 
larger  moments  in  disagreement  with  the  W3P  (  to-3  falloff  beyond  the 
comer  frequency)  model.  Also  note  the  explosion  spectral  amplitudes  gen¬ 
erally  fall  below  the  Sharpe  curve. 


the  explosion  data  points  have  less  30-Hz  energy  than  the  Sharpe 
model  predicts  and  fall  in  the  region  of  the  earthquake  population. 

Discussion 

Evermlcn  el  al.,  [  1986]  argue  that  the  Sharpe  model  is  a  conser¬ 
vative  or  minimum  estimate  for  the  high-frequency  energy  radi¬ 
ated  by  explosions.  They  expect  all  explosion  data  to  plot  on  or 
above  this  curve.  The  reason  for  die  failure  of  all  the  explosion 
data  to  plot  on  or  above  the  Sharpe  curve  may  be  due  to  loss  of 
high  frequency  energy  in  the  transmission  from  source  to  receiver, 
in  the  source  spectrum  itself,  or  some  combination  of  both.  The 
transmission  characteristics  of  the  medium  at  these  high  frequen¬ 
cies  arc  not  known  well  enough  to  determine  the  magnitude  of 
this  effect.  Ideally  a  series  of  small  nuclear  explosions  in  hard 
rock  recorded  in  hard  rock  would  be  the  best  comparison  with  the 
cardtquakc  data  presented  here.  Since  these  data  are  not  available, 
we  will  be  confined  to  making  some  general  remarks  about  explo¬ 
sion  models. 

The  comer  frequency  and  falloff  at  high  frequencies  of  explo¬ 
sion  spectra  are  critical  to  the  Evernden  cl  al.  discrimination  pro¬ 
posal.  The  question  of  whether  the  step  function  in  pressure 
Sharpe  model  represents  a  minimum  estimate  of  the  high- 
frequetkx  energy  radiated  by  explosions  must  be  examined  in 
light  of  the  data  presented.  It  is  not  clear,  particularly  for  tamped 
explosions  at  shallow  depths  and  high  frequencies,  that  the  pres¬ 


sure  at  the  clastic  radius  is  a  step  function.  The  nonlinear 
processes  such  as  crushing  and  plastic  deformation  might  be 
expected  to  preferentially  attenuate  the  high-frequency  com¬ 
ponents  of  the  pressure  function.  Thus  it  seems  possible  tha'  a 
tamped  explosion  could  have  a  spectral  falloff  steeper  than  m  2 
and  have  less  energy  at  high  frequencies  than  the  Sharpe  model 
predicts.  Much  more  work  needs  to  be  done  in  determining  llie 
exact  high-frequency  behavior  of  both  tamped  and  decoupled 
explosions. 

Figure  1  implies  two  important  problems  for  high-frequency 
discrimination  regardless  of  the  actual  attenuation.  First,  the  high 
stress  drop  earthquakes  such  as  some  of  the  Oaxaca  aftershocks 
fall  close  to  or  above  the  Sharpe  curve.  In  fact  both  the  W2P  and 
W3P  models  predict  Uiat  high  stress  drop  earthquakes  can  radiate 
more  high-frequency  energy  dian  the  Sharpe  model.  Thus  the 
question  for  discrimination  is  not  only  whether  the  W2P  or  Vv'3!‘ 
better  describe  the  data  but  also  whether  high  stress  drop  earth 
quakes  are  expected  in  the  regions  of  interest  and  whcdicr  they 
can  be  discriminated  by  methods  other  than  high-frequency  spec¬ 
tral  amplitude. 

Second,  it  may  be  possible  to  defeat  any  high-frequency 
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Fig.  2.  Observed  and  theoretical  velocity  seismograms  passed  through  a 
narrow-band  filter  centered  at  30  Hz  (20-40  Hz).  The  Mammoth  trace 
represents  a  complex  earthquake.  The  two  Anza  traces  represent  a  single 
earthquake,  showing  simple  explosionlike  earthquake  signal  at  one  station 
(FKD)  and  a  more  typical  simple  seismogram  with  a  large  S  wave  pulse  at 
another  station  (PFO).  The  trace  for  the  Sharpe  explosion  model  has  all 
the  higher-  frequency  energy  arriving  in  the  impulsive  beginning,  whereas 
for  the  W3P  model  the  amplitude  of  the  stopping  phase  is  about  3  limes 
larger  than  the  starting  phase.  For  the  constant  rupture  velocity  and  stress 
drop  W3P  model,  the  30-Hz  spectrum  docs  not  increase  with  moment 
beyond  1020  dyn  cm  because  the  pulses  seen  through  a  narrow-band  filter 
only  spread  out  in  time,  but  do  not  change  amplitude. 


discrimination  technique  by  detonating  the  shot  in  a  medium  with 
high  attenuation  such  as  is  expected  at  Pahuie  Mesa.  Certainly  any 
technique  based  on  simple  spectral  amplitude  at  30  Hz  would  fail 
for  some  of  the  recordings  of  the  explosion  data  presented  here, 
even  for  events  recorded  as  close  as  5  km.  Scattering  and  attenua¬ 
tion  can  be  expected  to  have  a  larger  effect  on  the  source  spectrum 
at  these  high  frequencies. 

The  data  in  Figure  lc  show  that  large  earthquakes  radiate  more 
30-Hz  energy  than  small  ones,  in  contrast  with  a  constant  stress 
drop  W3P  model.  This  is  obvious  from  a  cursory  observation  of 
the  long-duration,  higher-amplitude  P  waves  from  larger  events. 
As  a  simple  illustration  of  this  point,  we  present  in  Figure  2, 
narrow-band  filtered  velocity  seismograms  (bandwidth  20-40  Hz) 
from  throe  earthquake  recordings  along  with  filtered  seismograms 
of  the  Sharpe  explosion  model  and  the  W3P  [Archambeau,  1968, 
1972)  earthquake  model.  These  records  clearly  indicate  tire 
greater  amount  of  30-Hz  energy  radiated  from  larger  events  due  to 
both  larger  peak  amplitudes  and  longer,  more  complex  time  func¬ 
tions.  In  contrast  the  W3P  model  in  Figure  2  does  not  radiate  any 
more  energy  at  30  Hz  as  the  moment  increases,  since  increasing 
the  source  size  only  causes  the  starting  and  stopping  phases  of  the 
W3P  model  to  separate  further  in  lime,  but  does  not  change  their 
amplitude. 

The  difference  in  the  time  domain  of  the  theoretical  models 
suggests  that  the  waveforms  themselves  be  considered  in  any 
high-frequency  discrimination  technique.  The  high-frequency 
energy  radiated  by  earthquakes  is  generated  by  rupture  complexity 
which  can  be  modeled  by  accelerations  and  decelerations  of  the 
rupture  front  [e  g.,  Bernard  and  Madariaga,  1984],  The  W3P 
model,  as  indicated  in  Figure  2,  radiates  high-frequency  energy  in 
a  starting  and  stopping  phase  with  the  latter  having  about  3  times 
higher  amplitude  than  the  former.  Explosions,  in  contrast,  arc 
expected  to  radiate  high-frequency  energy  in  an  initial  impulse,  as 
indicated  by  the  Sharpe  model  in  Figure  2.  In  some  cases,  earth¬ 


quake  records  show  an  impulsive  "explosionlike"  record  at  one  or 
two  stations  of  an  array,  but  in  these  cases,  the  other  stations  of 
the  array  usually  show  a  less  impulsive,  more  "earlhquakelike" 
signals  with  large  S  waves  (Figure  2).  There  may  be  a  problem, 
however,  in  applying  such  time  domain  or  phase  characteristic 
approach  to  high  stress  drop  earthquakes.  High  stress  drop  earth¬ 
quakes  have  high  comer  frequencies  and  smaller  source  dimen¬ 
sions  and  thus  shorter  rupture  times.  A  magnitude  3  or  4  earth¬ 
quake  with  a  stress  drop  of  100  bars  or  more  might  have  a  rupture 
duration  of  less  than  a  tenth  of  a  second,  and  might  make  discrim¬ 
ination  based  on  either  waveform  in  the  time  domain  or  spectral 
amplitude  difficult. 


Conclusions 

Observations  of  high-frequency  P  wave  spectra  from  earth¬ 
quakes,  normalized  to  a  distance  of  10  km,  show  an  increase  in 
spectral  amplitude  with  moment  approximately  in  agreement  with 
an  co~2  falloff  model,  and  for  large  moments  in  disagreement  with 
an  to"3  falloff  model.  The  explosion  data  at  a  frequency  of  30  Hz 
generally  fall  below  the  Sharpe  curve,  show  much  scatter,  and  do 
not  all  discriminate  from  the  earthquake  data  on  the  basis  of  spec¬ 
tral  amplitude.  How  much  of  this  effect  is  due  to  attenuation 
requires  further  investigation.  This  suggests  that  any  discrimina¬ 
tion  technique  using  high-frequency  P  waves  should  also  take  into 
account  phase  information  (e.g.,  time  domain  character),  not  just 
spectral  level.  High  stress  drop  earthquakes  can  radiate  more 
high-frequency  energy  than  the  Sharpe  model  and  might  be 
expected  to  make  any  technique  based  on  time  domain  character 
of  P  waves  difficult  as  well. 
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Abstract.  Broadband  receiver  functions  determined  from  teleseismic 
P- waveforms  recorded  at  two  seismic  stations  in  eastern  Kazakh, 
U  S  S  R.,  are  inverted  for  the  vertical  velocity  structure  beneath  the  sta¬ 
tions,  The  detailed  broadband  receiver  functions  are  obtained  by  stacking 
source-equalized  radial  components  of  teleseismic  P-waveforms.  A  time- 
domain  inversion  of  the  radial  receiver  function  is  used  to  determine  the 
structure  assuming  a  crustal  model  parameterized  by  flat-lying,  homo¬ 
geneous  layers.  The  general  features  of  the  inversion  results  are:  a  com¬ 
plex  shallow  crust,  velocities  !»"'  dian  6  km/s  in  the  upper  crust,  a  high 
velocity  (6.9-7.5  km/s)  lower  crust,  and  a  Moho  dial  varies  between  47 
and  57  km  depth.  These  results  compare  favorably  with  a  composite  velo¬ 
city  model  from  Deep  Seismic  Sounding  data,  but  show  lower  velocities 
in  the  upper  crust,  and  higher  velocities  in  the  lower  crust.  The  results 
indicate  that  the  crust  in  this  region  is  relatively  uniform,  however  the 
nature  and  depth  of  the  Moho  changes  significantly  across  the  region,  The 
receiver  function  inversion  structure  for  the  Kazakh  sites  is  similar  to  a 
published  receiver  function  structure  for  a  site  on  the  Canadian  Shield. 


Introduction 

Since  July,  1986  three  seismic  stations  have  been  operating  in  the  vicin¬ 
ity  of  the  Soviet  nuclear  test  site  in  eastern  Kazakh,  U.S.S.R,,  under  a  joint 
agreement  between  the  Natural  Resources  Defense  Council  (NR  DC)  and 
the  Soviet  Academy  of  Science  [Berger  et  al,  1987).  Data  from  July.  1986 
through  February,  1987  were  from  seismometers  placed  in  surface  vaults; 
since  March,  1987,  data  also  have  been  recorded  from  seismometers 
located  in  100  meter  deep  boreholes  beneath  the  surface  vaults.  All  data 
are  digitally  recorded.  Figure  1  shows  the  location  of  the  three  stations 
(KKL,  BAY,  KSU),  the  locations  of  teleseismic  earthquakes  used  in  this 
study,  and  the  Deep  Seismic  Sounding  (DSS)  line  discussed  by  Antonenko 
[1984], 

Leith  [written  communications,  1987]  has  compiled  geophysical  and 
geological  data  for  the  NRDC  sites,  from  the  Soviet  literature.  The 
seismograph  sites  lie  within  the  Kazakh  fold  system,  a  belt  of  folded  an .1 
faulted  Paleozoic  rocks.  The  region  was  marked  by  extensive  igneou 


Fig.  1 ,  The  lower  part  of  the  figure  shows  the  location  of  the  NRDC  seismic  stations,  Dcgelen  Mountain  on  the 
Soviet  test  site,  and  the  DSS  line  discussed  by  Antonenko  (1984).  The  upper  part  of  the  figure  shows  the  location  of 
the  study  area  in  central  Asia,  and  the  location  of  the  teleseismic  events  studied:  (■ )  denote  events  recorded  at 
KKL,  (•)  denote  events  recorded  at  BAY. 
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activity  during  the  late  stages  of  the  Hercynian  Orogeny  durinc  ■V 
Permo-Triassic  Period.  The  three  seismograph  sitrs  arc  located  or  s  r j 
granitic  intrusions  of  Paleozoic  to  early  Mesozoic  age  The  grarres  c 
intruded  into  deformed  Paleozoic  sedimentary  rocks.  Subsequent  tectui  . 
activity  has  resulted  in  regional  crustal  uplift  and  thickening  Lrub  ug 
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gests  that  the  geological  and  geophysical  characteristics  of  Degelen 
Mountain  on  the  Soviet  Test  Site  are  similar  to  those  of  the  Karkaralinsk 
and  Bayanaui  intrusions.  The  plutons  and  Degelen  Mountain  are  com¬ 
posed  of  alaskite,  an  alkali  feldspar  granite  of  slightly  lower  density  and 
seismic  velocity  than  typical  granite.  Based  on  Soviet  deep  seismic 
sounding  lines  in  the  area,  the  Pn-velocity  appears  to  vary  from  8.0  to  8.3 
km/s  and  the  crustal  thickness  vanes  from  50  km  in  the  vicinity  of  Kar¬ 
karalinsk  and  Bayanaui  to  possibly  as  little  as  40  km  at  Karasu.  Soviet 
geophysical  results  indicate  that  at  least  the  Karkaralinsk  intrusion  is  a  lac¬ 
colith  narrowing  from  an  area  of  about  450  km2  at  the  surface  to  about 
100  km2  at  a  depth  of  10  km. 

In  this  report,  we  model  receiver  functions  obtained  from  leleseismic 
P-waveforms  recorded  at  two  sites,  Karkaralinsk  and  Bayanaui,  to  deter¬ 
mine  the  detailed  structure  of  the  crust  and  upper  mantle  directly  beneath 
the  recording  site.  Insufficient  data  has  been  recorded  from  the  Karasu 
site  for  this  type  of  analysis.  The  broadband  receiver  function  method  is 
most  sensitive  to  the  details  of  the  vertical  shear  velocity  distribution  in 
the  ciusi.  a  difficult  parameter  to  constrain  using  other  methods.  We  esti¬ 
mate  gross  lateral  variations  in  structure  by  examining  the  variation  in 
vertical  velocity  distribution  azimuthally  around  the  station,  and  by  exa¬ 
mining  the  relative  magnitude  of  the  transverse  and  radial  arrivals. 
Finally,  assuming  a  constant  Poisson's  ratio  we  compare  our  results  with 
P-velocity  profiles  from  a  nearby  DSS  line. 

Data  and  Analysis 

Wc  have  used  the  deconvolution  procedure  proposed  by  Langston 
(1979|  to  isolate  the  response  of  the  structure  beneath  the  seismograph 
(the  receiver  function)  from  the  source  and  path  effects,  all  of  which 
interact  to  form  the  observed  seismograms  recorded  at  leleseismic  dis¬ 
tances.  For  teleseisms  recorded  at  30°  to  90°  epicentre!  distance,  the  initial 
energy  arriving  from  the  source  region  consists  of  plane  P-wave  arrivals. 
Since  the  earth  structure  beneath  the  station  will  produce  P  lo  S  conver¬ 
sions  the  horizontal  component  of  ground  motion  will  in  general  be  quite 
different  from  the  vertical  component.  Langston  [1979]  has  shown  that 
the  receiver  function  can  be  obtained  by  deconvolving  the  vertical  com¬ 
ponent  from  the  radial.  The  resulting  receiver  functions  from  events 
clustered  both  in  distance  and  back  azimuth  are  then  slacked  m  the  lime 
domain  to  obtain  a  single  estimate  of  the  horizontal  receiver  function. 
Stacking  the  deconvolved  traces  from  individual  events  averages  small 
differences  in  the  estimated  mean  receiver  functions  and  improves  the 
signal-to-noisc  ratios. 

Figure  2  shows  seismograms  recorded  at  KKL  from  shallow  and  deep 
focus  events  in  Japan.  The  upper  three  traces  are  the  vertical,  radial  and 
transverse  components,  respectively.  The  recorded  seismograms  show  the 
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Fig  2  Deep  and  shallow  focas  events  recorded  at  KKL  The  upper  three 
traces  are  the  vertical,  radial,  and  transverse  components,  respectively 
The  lower  two  traces  arc  the  deconvolved  radial  and  transverse  receiver 
functions.  The  horizontal  components  arc  plotted  at  the  same  scale 


differences  in  complexity  typically  seen  between  deep  and  shallow  focus 
events.  The  bottom  two  traces  are  the  deconvolved  radial  and  transverse 
components.  The  deconvolved  radial  components  for  both  events  show 
very  similar  features,  although  these  two  traces  exhibit  the  largest  differ¬ 
ences  among  individual  seismograms  used  in  any  of  the  three  stacks 
reported  in  this  paper.  The  seismograms  available  were  sufficient  to  pro¬ 
duce  three  stacked  receiver  functions:  Two  for  separate  azimuths  for  KKL 
and  one  for  BAY.  The  stack  for  KKL- Japan  includes  three  events;  the 
stack  for  KKL-Sumatra  has  4  events;  and  die  BAY-Japan  stack  has  6 
events.  The  stacked  radial  receiver  functions  lend  to  decay  rapidly  in 
amplitude,  indicating  a  relatively  homogeneous  crust,  "be  transverse 
components  are  generally  small  and  have  nearly  constant  plnude.  indi¬ 
cating  that  lateral  variations  in  the  structure  are  small. 

The  stacked  receiver  functions  were  inverted  using  a  standard,  linear¬ 
ized  inversion  procedure  [Owens  et  al.  1984]  with  the  modified  fast  partial 
derivatives  method  of  Randall  [1988],  The  earth  structure  is  parameter¬ 
ized  by  a  series  of  plane  horizontal  layers,  and  an  iterative  procedure  is 
used  to  minimize  the  residual  difference  between  the  radial  receiver  func¬ 
tion  and  the  synthetic  predicted  by  the  structure.  Our  starting  model  for 
the  KKL  and  BAY  inversions  is  from  the  structure  compiled  by  Leith 
[written  communications,  1987]  from  the  Soviet  DSS  literature.  The 
inversion  model  is  parameterized  with  1 .0  km  thick  layers  in  the  shallow- 
crust  (upper  5  km),  and  2.5  km  thick  layers  to  a  depth  of  57.5  Icm.  We 
invert  for  the  vertical  shear  wave  velocity  structure  for  each  back  azimuth 
for  which  we  had  sufficient  data:  the  compressional  velocities  arc  oblaincd 
by  assuming  a  Poisson’s  ratio  of  0.25.  and  die  densities  are  obtained  using 
the  relationship  r  =  0.32  Vp  +  0.77  (Owens  et  al .  1984)  The  vanauon  in 
structure  with  azimuth  allows  a  qualitative  esumaie  of  the  lateral  vanauon 
in  structure  about  the  station.  The  results  bom  the  thin-layer  inversions 
were  examined  to  identify  the  depths  of  the  significant  seismic  impedance 
contrasts,  then  the  models  were  reparameterized  with  thicker  layers  and 
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Fig.  3.  Companion  of  the  starling  model  and  final  inversion  results  The 
inversion  results  are  in  terms  of  the  shear  wave  structure,  however  for 
comparison  with  the  DSS  results,  the  models  are  presented  as  compres¬ 
sional  wave  velocity  assuming  a  Poisson's  rauo  of  0  25. 
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Idblc  i 

1  injl  Velocity  Models 

Leith 
h  a 

BAY-Japan 
h  a 

KKL-Japan 
h  a 

KKL-Sumatra 
h  a 

1.0  5.10 

1.0  5.32 

1.0  5.19 

1.0  5.14 

4.0  5.60 

3.0  5.96 

4.0  5.85 

4.0  5.91 

7.0  6.05 

6.0  5.66 

2.5  5.77 

7.5  5.51 

12.0  6.35 

5.0  6.01 

2.5  5.57 

2.5  6.03 

12.0  6.60 

10.0  6.29 

2.5  5.90 

12.5  6.51 

7.0  6.90 

5.0  6.65 

5.0  6.19 

12.5  7.08 

7.0  7  00 

15.0  6.89 

10  0  6  52 

2.5  7.29 

~>  8.35 

2.5  7.59 

15.0  7.04 

15.0  7.49 

~  8.33 

10.0  7.47 

5.0  7.80 
°°8  10 

•o  8.10 

h  =  layer  thickness  in  km,  a  =  P-velocity  in  km/s 


inverted  again.  All  the  inversions  used  in  this  study  converged  and  final 
results  included  all  singular  values;  hence  the  model  resolution  is  theoreti¬ 
cally  perfect  We  attempted  to  determine  final  models  with  the  minimum 
number  of  parameters  to  obtain  a  satisfactory  fit  to  the  data. 

Results 

The  final  inversion  velocity  models  are  shown  in  Figure  3  and  tabulated 
in  Table  1.  The  models  are  shown  in  terms  of  the  P-wave  velocity  struc¬ 
ture,  assuming  a  Poisson's  ratio  of  0.25.  Leith's  summary  of  the  DSS 
results,  which  was  used  as  a  starting  model  for  the  inversions,  is  shown  as 
the  solid  line.  Two  models  are  shown  for  KKL,  one  for  data  from  Japan 
earthquakes  and  one  from  Sumatra  earthquakes.  The  2a  error  bounds 
varied  with  depth  but  averaged  about  0.1  km/s,  however  the  differences 
between  the  two  KKL  models  may  be  more  representative  of  the  model 
variance.  A  single  BAY  structure  is  shown  for  data  from  Japan  earth¬ 
quakes. 

The  data  from  the  three  stacks  were  inverted  separately  with  a  model 
parameterized  by  thin  layers  and  using  the  Leith  starting  model.  The 
resulting  inversion  models  have  a  complex  shallow  crust  with  thin  high- 
velocity  layers.  The  upper  few  kilometers  of  the  structures  are  compli¬ 
cated  by  the  laccolith  intrusions  and  thus  may  be  expected  to  deviate  from 
the  plane  layered  structure  assumed.  Also,  pronounced  side  lobes  about 
die  main  arrival  (sec  Figure  2)  result  from  deconvolution  of  the  band- 
limited  data  (the  synthetic  seismograms  were  treated  in  the  same  manner), 
and  some  of  the  complexity  in  the  shallow  crust  results  from  the  inversion 
trying  to  fit  this  artifact.  To  assess  the  effect  of  this  in  the  deeper  structure 
we  have  recomputed  the  synthetics  with  a  smoother  shallow  crust  and 
thicker  layer  parameters,  as  we  described  earlier.  The  subsequent  synthet¬ 
ics  still  fit  the  lower  frequency  component  of  the  data,  indicating  that  the 
overall  features  of  the  structure  are  not  biased  by  the  complexities  in  the 
shallow  crust. 

The  overall  features  of  the  two  KKL  models  for  data  from  two  azimuths 
about  75°  apart  arc  the  same,  implying  no  pronounced  systematic 
nonplanar  structure.  This  is  supported  by  the  small  signals  in  the  decon¬ 
volved  tangential  components  (Fig.  2).  Compared  with  the  summary  of 
the  DSS  data  from  Leith,  velocities  in  the  upper  crust  (5-15  km  depth)  are 
slower,  less  than  6  km/s.  The  velocity  model  from  inversion  of  the  KKL 
receiver  function  is  faster  in  the  lower  crust  compared  to  the  DSS  model, 
with  velocities  reaching  7  km/s  at  depths  as  shallow  as  27  5  km  as 
opposed  to  36  km  depth  in  the  DSS  model  The  most  significant  differ¬ 
ence  is  in  the  nature  of  the  Moho  The  DSS  model  has  a  very  pronounced 
Moho  at  50  km  depth,  whereas  the  receiver  function  model;;  have  a  much 
less  distinct  Moho  at  52-57  km  The  upper  two  traces  of  Figure  4  show 
the  fits  of  the  synthetic  radial  receiver  functions  calculated  for  the  two 
KKL  models  of  Figure  3,  ovcrplotlcd  on  the  radial  receiver  function 
stacks  The  fits  for  the  thin  layer  inversions  were  much  better,  however 
we  did  not  feel  justified  in  fitting  the  data  to  that  extent,  and  feel  the 
coarser  models  shown  in  Figure  1  are  more  realisuc. 


Fig.  4.  Fit  of  synthetic  receiver  functions  (dotted  lines)  and  the  observed 
(solid  lines)  stacked  radial  receiver  functions.  Top  traces  for  KKI 
Sumatra,  middle  for  KJCL-Japan,  and  bottom  for  BAY-Japan. 

The  result  for  the  BAY-Japan  inversion  is  also  shown  in  Figure  7  and 
tabulated  in  Table  1,  and  the  fit  of  the  synthetic  to  the  stack  is  show  n  in 
Figure  4.  In  general  the  features  of  the  BAY  model  are  similar  to  those 
for  KKL:  a  complex  shallow  crust,  an  upper  crustal  low  velocity  zone  and 
a  mid-crustal  interface  near  15  km;  however  the  lower  crustal  velocitie: 
are  generally  slower  and  the  Moho  is  both  shallower  and  more  distinct 
than  observed  at  KKL.  In  fact,  below  10  km  the  BAY  velocity  structure  is 
quite  close  to  die  DSS  velocity  structure  except  the  Moho  is  shallower  (45 
-  47.5  km  depth)  than  in  the  DSS  model  (50  km).  Although  we  have  only 
inverted  BAY  data  from  one  azimuth,  the  small  signal  observed  on  the 
deconvolved  tangential  component  relative  to  the  radial  implies  that  the 
Structure  does  not  depart  systematically  from  the  assumed  horizontal 
layers. 

The  largest  differences  among  all  three  inversion  models  arc  in  die 
depth  and  nature  of  the  Moho.  The  Moho  from  the  KKL -Japan  inversion 
is  less  pronounced  than  that  from  the  JCKL-Sumalra  inversion.  The  rela¬ 
tively  strong  arrival  in  die  KKL-Sumatra  stack  at  7  sec  is  predominantly 
the  Moho  P-lo-S  conversion  and  is  the  controlling  factor  for  the  depth  ar.-’ 
impedance  contrast  of  that  interface.  The  arrival  is  weak  or  missing  in  the 
KKL -Japan  stack.  This  azimuthal  variation  may  reflect  variations  in  the 
amplitude  of  a  converted  phase  from  a  dipping  layer  or  the  interference 
between  multiple  arrivals  which  can  be  easily  destroyed  by  lateral  varia¬ 
tions  in  velocity  or  layer  thickness.  The  lack  of  corresponding  systematic 
arrivals  on  die  transverse  component  favors  die  interference  explanation; 
also  the  Soviet  DSS  studies  in  the  region  indicate  significant  lateral  varia¬ 
tions  in  the  Moho.  The  prominent  Ps  phase  in  the  BAY-Japan  stack  al  5.5 
seconds  controls  the  shallower  depth  and  larger  contrast  of  the  Moho  on 
die  BAY  model.  Zandl  and  Owens  (1986)  discussed  a  similar  companson 
between  refraction  and  receiver  function  derived  velocity  structure* 
Cumberland  Plateau,  eastern  Tennessee.  The  more  distinct  nature  t  ■' ' 
refraction  Moho  compared  with  the  receiver  function  Moho  obseivu*  r. 
both  studies  most  likely  results  from  differences  in  the  frequency  a  .J 
phase  velocities  of  die  arrivals  utilized  by  the  two  techniques 

Summary  and  Conclusions 

Tclescismic  receiver  functions  from  broadband  waveforms  recorded  ir, 
eastern  Kazakh  have  been  inverted  for  the  crustal  vclocuv  structure  xml 
compared  with  a  summary  of  DSS  data  compiled  by  Leith  The  results  an 
summarized  in  terms  of  P-wave  velocity  models  in  Table  !  and  plotted  in 
Figure  3.  Our  results  are  in  genera)  agreement  with  the  DSS  1  a  ■ 
although  our  models  indicate  the  presence  of  a  high  vcloc’ty  (>  '  7  •  ■  i  ' 
layer  below  30  km  which  is  not  seen  in  the  DSS  results  This  t.-pr.-  i 
either  real  differences  in  velocity  or  variations  in  Poisson's  •aiu  iV 
inversion  structures  from  data  recorded  at  the  two  sues,  wh.Oi  ,m-  11  1 
imalely  150  km  apart,  are  similar  tndinuing  that  the  crust  m  thn  - 
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Vp  (km/aac) 


*  *  » 


Fig.  5.  Comparison  of  the  Kazakh  receiver  function  inversion  models 
with  the  receiver  function  inversion  model  of  Owens  et  al  (1987)  for  the 
RSON  site  on  the  Canadian  Shield. 


may  be  relatively  uniform.  However,  the  nature  of  the  Moho  is  different 
between  the  two  sites;  at  Bayanaul  it  is  shallower  (47  km  depth)  and  much 
more  distinct  than  at  Karkaralinsk  (52-57  km  depth).  The  shallow  crust  is 
complex,  showing  thin  high  velocity  layers  above  an  upper  crustal  low 
velocity  zone.  These  shallow  features  may  not  be  regional,  but  may  be 
largely  confined  to  the  laccoliths  in  which  the  seismic  stations  are  sited. 

In  Figure  5  the  Bayanaul  and  Karkaralinsk  receiver  function  models  arc 
compared  with  a  smoothed  receiver  function  model  for  Red  Lake,  Ontario 
(RSON)  on  the  Canadian  Shield  [Owens  ct  al,  1987).  Studies  of  data  from 
RSON  have  often  been  used  as  analogs  for  the  Kazakh  region  of  the 
Soviet  Union.  The  shallow  crust  at  RSON  has  higher  velocities  than  our 
Kazakh  models,  however  this  may  be  from  complexities  introduced  from 
the  Kazakh  stations  being  sited  in  the  laccolith  intrusions,  whereas  RSON 


is  located  on  a  more  uniform  upper  crust  Also,  the  Kazakh  sites  are 
located  in  aiaskite  granitic  intrusions  which  have  lower  velocity  than  the 
metasedimentary-  gneiss  terrain  where  RSON  is  sited.  Between  10  and  40 
km  depth  the  Kazakh  and  RSON  models  are  nearly  identical.  The  RSON 
model  has  a  somewhat  lower  velocity  in  the  lower  crust  (26  to  36  km 
depth),  but  has  a  shallower  Moho  (40  to  45  km  depth)  than  observed  for 
Kazakh. 
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10210  Campus  Point  Drive 
San  Diego,  CA  92121  (2  copies) 


J.  Barker 

Department  of  Geological  Sciences 
State  University  of  New  York 
at  Binghamton 
Vestal,  NY  13901 

Dr.  T.J.  Bennett 
S -CUBED 

A  Division  of  Maxwell  Laboratory 
1 1800  Sunrise  Valley  Drive,  Suite  1212 
Rcston,  VA  22091 

Mr.  William  J.  Best 
907  Westwood  Drive 
Vienna,  VA  22180 


Dr.  N.  Biswas 
Geophysical  Institute 
University  of  Alaska 
Fairbanks,  AK  99701 


Dr.  G.A.  Bollinger 
Department  of  Geological  Sciences 
Virginia  Polytechnical  Institute 
21044  Derring  Hall 
Blacksburg,  VA  24061 

Dr.  Stephen  Bratt 
Science  Applications  Int'l  Corp. 
10210  Campus  Point  Drive 
San  Diego,  CA  92121 


Michael  Browne 
Teledyne  Geotech 
3401  Shiloh  Road 
Garland,  TX  75041 


Mr.  Roy  Burger 
1221  Serry  Road 
Schenectady,  NY  12309 


Dr.  Robert  Burridge 
Schlumberger-Doll  Research  Center 
Old  Quarry  Road 
Ridgefield,  CT  06877 


Dr.  Jerry  Carter 
Rondout  Associates 
P.O.  Box  224 
Stone  Ridge,  NY  12484 


Dr.  W.  Winston  Chan 
Teledyne  Geotech 
314  Montgomery  Street 
Alexandria,  VA  22314-1581 


Dr.  Theodore  Cherry 
Science  Horizons,  Inc. 

710  Encinitas  Blvd.,  Suite  200 
Encinitas,  CA  92024  (2  copies) 


Prof.  Jon  F.  Claerbout 
Department  of  Geophysics 
Stanford  University 
Stanford,  CA  94305 


Prof.  Robert  W.  Clayton 

Scismological  Laboratory 

Division  of  Geological  &  Planetary  Sciences 

California  Institute  of  Technology 

Pasadena,  CA  91125 


Prof.  P.  A.  Dahlen 

Geological  and  Geophysical  Sciences 

Princeton  University 

Princeton,  NJ  08544-0636 


Prof.  Anton  W.  Dainty 
Earth  Resources  Lab 
Massachusetts  Institute  of  Technology 
,  42  Carleton  Street 

Cambridge,  MA  02142 

Prof.  Adam  Dziewonski 
Hoffman  Laboratory 
Harvard  University 
20  Oxford  St 
Cambridge,  MA  02138 

Prof.  John  Ebel 

Department  of  Geology  &  Geophysics 
Boston  College 
Chestnut  Hill,  MA  02167 


Eric  Fielding 
SNEE  Hall 
INSTOC 

Cornell  University 
Ithaca,  NY  14853 

Prof.  Donald  Forsyth 
Department  of  Geological  Sciences 
Brown  University 
Providence,  RI  02912 


Prof.  Art  Frankel 
Mail  Stop  922 
Geological  Survey 
790  National  Center 
Reston,  VA  22092 

Dr.  Anthony  Gangi 
Texas  A&M  University 
Department  of  Geophysics 
College  Station,  TX  77843 


Dr.  Freeman  Gilbert 

Inst,  of  Geophysics  &  Planetary  Physics 
University  of  California,  San  Diego 
P.O.  Box  109 
La  Jolla,  C A  92037 

Mr.  Edward  Giller 
Pacific  Sierra  Research  Corp. 

1401  Wilson  Boulevard 
Arlington,  VA  22209 


Dr.  Jeffrey  W.  Given 
Sierra  Geophysics 
11255  Kirkland  Way 
Kirkland,  WA  98033 


Prof.  Stephen  Grand 
University  of  Texas  at  Austin 
Department  of  Geological  Sciences 
Austin,  TX  78713-7909 


Prof.  Roy  Greenfield 
Geosciences  Department 
403  Deike  Building 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

Dan  N.  Hagedom 
Battelle 

Pacific  Northwest  Laboratories 
Battelle  Boulevard 
Richland,  WA  99352 

Kevin  Hutchenson 
Department  of  Earth  Sciences 
St.  Louis  University 
3507  Laclede 
St.  Louis,  MO  63103 

Prof.  Thomas  H.  Jordan 
Department  of  Earth,  Atmospheric 
and  Planetary  Sciences 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 

Robert  C.  Kemerait 
ENSCO,  Inc. 

445  Pineda  Court 
Melbourne,  FL  32940 


William  Kikendall 
Teledyne  Geotech 
3401  Shiloh  Road 
Garland,  TX  75041 


Prof.  Leon  Knopoff 

University  of  California 

Institute  of  Geophysics  &  Planetary  Physics 

Los  Angeles,  CA  90024 


Prof.  L.  Timothy  Long 
School  of  Geophysical  Sciences 
Georgia  Institute  of  Technology 
Atlanta,  GA  30332 
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Prof.  Art  McGarr 
Mail  Stop  977 
Geological  Survey 
345  Middlefield  Rd. 
Menlo  Park,  CA  94025 

Dr.  George  Mellman 
Sierra  Geophysics 
1 1255  Kirkland  Way 
Kirkland,  WA  98033 


Prof.  John  Nabelek 
College  of  Oceanography 
Oregon  State  University 
Corvallis,  OR  97331 


Prof.  Geza  Nagy 

University  of  California,  San  Diego 
Department  of  Ames,  M.S.  B-010 
La  Jolla,  CA  92093 


Prof.  Amos  Nur 
Department  of  Geophysics 
Stanford  University 
Stanford,  CA  94305 


Prof.  Jack  Oliver 
Department  of  Geology 
Cornell  University 
Ithaca,  NY  14850 


Prof.  Robert  Phinney 
Geological  &  Geophysical  Sciences 
Princeton  University 
Princeton,  NJ  08544-0636 


Dr.  Paul  Pomeroy 
Rondout  Associates 
P.O.  Box  224 
Stone  Ridge,  NY  12484 


Dr.  Jay  Pulli 
RADIX  System,  Inc. 

2  Taft  Cour*1,  Suite  203 
Rockville,  MD  20850 


Dr.  Norton  Rimer 
S-CUBED 

A  Division  of  Maxwell  Laboratory 
P.O.  Box  1620 
La  Jolla,  CA  92038-1620 


Prof.  Larry  J.  Ruff 
Department  of  Geological  Sciences 
1006  C.C.  Little  Building 
University  of  Michigan 
Ann  Arbor,  MI  48109-1063 

Dr.  Richard  Sailor 
TASC  Inc. 

55  Walkers  Brook  Drive 
Reading,  MA  01867 


Thomas  J.  Sereno,  Jr. 

Science  Application  Int'l  Corp. 
10210  Campus  Point  Drive 
San  Diego,  CA  92121 


John  Sherwin 
Teledyne  Geotech 
3401  Shiloh  Road 
Garland,  TX  75041 


Prof.  Robert  Smith 
Department  of  Geophysics 
University  of  Utah 
1400  East  2nd  South 
Salt  Lake  City,  UT  84112 

Prof.  S.  W.  Smith 
Geophysics  Program 
University  of  Washington 
Seattle,  WA  98195 


Dr.  Stewart  Smith 
IRIS  Inc. 

1616  North  Fort  Myer  Drive 
Suite  1440 

Arlington,  VA  22209 

Dr.  George  Sutton 
Rondout  Associates 
P.O.  Box  224 
Stone  Ridge,  NY  12484 


Prof.  L.  Sykes 

Lamont-Doherty  Geological  Observatory 
of  Columbia  University 
Palisades,  NY  10964 


Prof.  Pradeep  Talwani 
Department  of  Geological  Sciences 
University  of  South  Carolina 
Columbia,  SC  29208 


Prof.  Ta-liang  Teng 
Center  for  Earth  Sciences 
University  of  Southern  California 
University  Park 
Los  Angeles,  CA  90089-0741 

-  Dr.  R.B.  Tittmann 

Rockwell  International  Science  Center 
1049  Camino  Dos  Rios 
,  P.O.  Box  1085 

Thousand  Oaks,  CA  91360 

Dr.  Gregory  van  der  Vink 
IRIS,  Inc. 

1616  North  Fort  Myer  Drive 
Suite  1440 

Arlington,  VA  22209 

William  R.  Walter 
Seismological  Laboratory 
University  of  Nevada 
Reno,  NV  89557 


Dr.  Gregory  Wojcik 
Weidlinger  Associates 
4410  El  Camino  Real 
Suite  1 10 

Los  Altos,  CA  94022 

Prof.  John  H.  Woodhouse 
Hoffman  Laboratory 
Harvard  University 
20  Oxford  Street 
Cambridge,  MA  02138 

Dr.  Gregory  B.  Young 
ENSCO,  Inc. 

5400  Port  Royal  Road 
Springfield,  VA  22151-2388 


Dr.  Cliff  Frolich 
Institute  of  Geophysics 
8701  North  Mopac 
Austin,  TX  78759 
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GOVERNMENT 


Dr.  Ralph  Alewine  III 
DARPA/NMRO 
1400  Wilson  Boulevard 
Arlington,  VA  01731-5000 


Mr.  James  C.  Battis 
GL/LWH 

Hanscom  AFB,  MA  22209-2308 


Paul  Johnson 

ESS-4,  Mail  Stop  J979 

Los  Alamos  National  Laboratory 

Los  Alamos,  NM  87545 


Janet  Johnston 
GL/LWH 

Hanscom  AFB,  MA  01731-5000 


Dr.  Robert  Blandfoid 
DARPA/NMRO 
1400  Wilson  Boulevard 
Arlington,  VA  87185 


Eric  Chael 
Division  9241 
Sandia  Laboratory 
Albuquerque,  NM  01731-5000 


Dr.  John  J.  Cipar 
GL/LWH 

Hanscom  AFB,  MA  01731-5000 


Mr.  Jeff  Duncan 
Office  of  Congressman  Markey 
2133  Rayburn  House  Bldg. 
Washington,  D.C.  20515 


Dr.  Jack  Evemden 
USGS  -  Earthquake  Studies 
345  Middlefield  Road 
Menlo  Park,  CA  94025 


Dr.  Katharine  Kadinsky-Cade 
GL/LWH 

Hanscom  AFB,  MA  01731-5000 


Ms.  Ann  Kerr 
IGPP,  A-025 

Scripps  Institute  of  Oceanography 
University  of  California,  San  Diego 
La  Jolla,  CA  92093 

Dr.  Max  Koontz 
US  Dept  of  Energy/DP  5 
Forrestal  Building 
1000  Independence  Avenue 
Washington,  DC  20585 

Dr.  W.H.K.  Lee 

Office  of  Earthquakes,  Volcanoes, 
&  Engineering 
345  Middlefield  Road 
Menlo  Park,  CA  94025 

Dr.  William  Leith 
U.S.  Geological  Survey 
Mail  Stop  928 
Reston,  VA  22092 


Art  Frankel 
USGS 

922  National  Center 
Reston,  VA  22092 


Dr.  T.  Hanks 
USGS 

Nat'l  Earthquake  Research  Center 
345  Middlefield  Road 
Menlo  Park,  CA  94025 

Dr.  James  Hannon 

Lawrence  Livermore  Nat'l  Laboratory 
P.O.  Box  808 
Livermore,  CA  94550 


Dr.  Richard  Lewis 

Director,  Earthquake  Engineering  &  Geophysics 
U.S.  Army  Corps  of  Engineers 
Box  631 

Vicksburg,  MS  39180 

James  F.  Lewkowicz 
GL/LWH 

Hanscom  AFB,  MA  01731-5000 


Mr.  Alfred  Lieberman 
ACDA/VI-OA’State  Department  Bldg 
Room  5726 
320  -  21st  Street,  NW 
Washington,  DC  20451 
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Stephen  Marigmu 
GL/LWH 

Manscom  AFB,  MA  01731-5000 


Dr.  Robert  Masse 
Box  25046,  Mail  Stop  967 
Denver  Federal  Center 
Denver,  CO  80225 


Art  McGarr 

U.S.  Geological  Survey,  MS-977 
345  Middlefield  Road 
Menlo  Park,  CA  94025 


Richard  Morrow 
A  CD  A/VI,  Room  5741 
320  21st  Street  N.W 
Washington,  DC  20451 


Dr.  Keith  K.  Nakanishi 
Lawrence  Livermore  National  Laboratory 
P.O.  Box  808,  L-205 
Livermore,  CA  94550 


Dr.  Carl  Newton 

Los  Alamos  National  Laboratory 

P.O.  Box  1663 

Mail  Stop  C335,  Group  ESS-3 
Los  Alamos,  NM  87545 

Dr.  Kenneth  H.  Olsen 

Ixis  Alamos  Scientific  Laboratory 

P.O.  Box  1663 

Mail  Stop  C335,  Group  ESS-3 
Los  Alamos,  NM  87545 

H  n/ord  J.  Patton 

Lawrence  Livermore  National  Laboratory 
P.O.  Box  808,  L-205 
Livermore,  CA  94550 


Dr.  Frank  F.  Pilotte 

HQ  AFTAOTT 

Patrick  AFB,  FL  32925-6001 


Katie  Poley 
CIA-OSWR/NED 
Washington,  DC  20505 


Mr.  Jack  Rachlin 
U.S.  Geological  Survey 
Geology,  Rm  3  Cl 36 
Mail  Stop  928  National  Center 
Reston,  VA  22092 

Dr.  Robert  Reinke 
WL/NTESG 

Kirtland  AFB,  NM  87117-6008 


Dr.  Byron  Ristvet 

HQ  DNA,  Nevada  Operations  Office 

Attn:  NVCG 

P.O.  Box  98539 

Las  Vegas,  NV  89193 

Dr.  George  Rothe 
HQ  AFTAC/TGR 
Patrick  AFB,  FL  32925-6001 


Dr.  Michael  Shore 
Defense  Nuclear  Agency/SPSS 
6801  Telegraph  Road 
Alexandria,  VA  22310 


Donald  L.  Springer 

Lawrence  Livermore  National  Laboratory 
P.O.  Box  808,  L-205 
Livermore ,  CA  94550 


Mr.  Chris  Paine 

Office  of  Senator  Kennedy,  SR  315 

United  States  Senate 
Washington,  DC  20510 

Colonel  Jerry  J.  Perrizo 
AFOSR/NP,  Building  410 
Bolling  AFB 

Washington,  DC  20332-6448 


Dr.  Lawrence  Turnbull 
OSWR/NED 

Central  Intelligence  Agency,  Room  5G48 
Washington,  DC  20505 


Dr.  Thomas  Weaver 
Los  Alamos  National  Laboratory 
P.O.  Box  1663,  Mail  Stop  C335 
Los  Alamos,  NM  87545 
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J.J.  Zucca 

Lawrence  Livermore  National  Laboratory 
Box  808 

Livermore,  CA  94550 


GL/SULL 
Research  Library 

Hanscom  AFB  ,  MA  01731-5000  (2  copies) 


Secretary  of  the  Air  Force  (SAFRD) 
Washington,  DC  20330 


Office  of  the  Secretary  Defense 
DDR&E 

Washington,  DC  20330 


HQ  DNA 

Attn:  Technical  Library 
Washington,  DC  20305 


DARPA/RMO/RETRIEVAL 
1400  Wilson  Boulevard 
Arlington,  VA  22209 


DARPA/RMO/Security  Office 
1400  Wilson  Boulevard 
Arlington,  VA  22209 


Geophysics  Laboratory 
Attn:  XO 

Hanscom  AFB,  MA  01731-5000 


Geophysics  Laboratory 
Attn:  LW 

Hanscom  AFB,  MA  01731-5000 


Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  VA  22314  (5  copies) 


Defense  Intelligence  Agency 
Directorate  for  Scientific  & 
Technical  Intelligence 
Washington,  DC  20301 


AFTAC/CA 

(STINFO) 

Patrick  AFB,  FL  32925-6001 


TACTEC 

Battelle  Memorial  Institute 
505  King  Avenue 

Columbus,  OH  43201  (Final  Report  Only) 


Mr.  Charles  L.  Taylor 
GL/LWH 

Hanscom  AFB,  MA  01731-5000 


DARPA/PM 

1400  Wilson  Boulevard 

Arlington,  VA  22209 
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Dr.  Ramon  Cabre,  SJ. 
Observatono  San  Calixto 
Casilla  5939 
La  Paz,  Bolivia 


w  Prof.  Hans-Peter  Harjes 
Institute  for  Geophysik 
Ruhr  University/Bochum 
»  P.O.  Box  102148 
4630  Bochum  1,  FRG 

Prof.  Eystein  Husebye 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjeller,  NORWAY 


Prof.  Brian  L.N.  Kennett 
Research  School  of  Earth  Sciences 
Institute  of  Advanced  Studies 
G.P.O.  Box  4 

Canberra  2601,  AUSTRALIA 

Dr.  Bernard  Massinon 

Societe  Radiomana 

27  rue  Claude  Bernard 

75005  Paris,  FRANCE  (2  Copies) 


Dr.  Pierre  Mecheler 
Societe  Radiomana 
27  rue  Claude  Bernard 
75005  Paris,  FRANCE 


Dr.  Svein  Mykkeltveit 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjeller,  NORWAY 


* 


CONTRACTORS  (Foreign) 
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FOREIGN  (Others) 


Dr.  Peter  Basham 

Earth  Physics  Branch 

Geological  Survey  of  Canada 

1  Observatory  Crescent 

Ottawa,  Ontario,  CANADA  K1A  0Y3 

Dr.  Eduard  Berg 
Institute  of  Geophysics 
University  of  Hawaii 
Honolulu,  HI  96822 


Dr.  Michel  Bouchon 
I.R.I.G.M.-B.P.  68 
38402  St.  Martin  D'Heres 
Cedex,  FRANCE 


Dr.  Hilmar  Bungum 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjeller,  NORWAY 


Dr.  Michel  Campillo 
Observatoire  de  Grenoble 
I.R.I.G.M.-B.P.  53 
38041  Grenoble,  FRANCE 


Dr.  Kin  Yip  Chun 
Geophysics  Division 
Physics  Department 
University  of  Toronto 
Ontario,  CANADA  M5S  1A7 

Dr.  Alan  Douglas 
Ministry  of  Defense 
Blacknest,  Brimpton 

Reading  RG7-4RS,  UNITED  KINGDOM 


Dr.  Roger  Hansen 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjeller,  NORWAY 


Dr.  Manfred  Henger 

Federal  Institute  for  Geosciences  &  Nat'l  Res. 
Postfach  510153 
D-3000  Hanover  51,  FRG 


Dr.  Fekadu  Kebede 
Seismological  Section 
Box  12019 

S-750  Uppsala,  SWEDEN 


Dr.  Tormod  Kvaema 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjeller,  NORWAY 


Dr.  Peter  Marshal 
Procurement  Executive 
Ministry  of  Defense 
Blacknest,  Brimpton 

Reading  FG7-4RS,  UNITED  KINGDOM 

Prof.  Ari  Ben-Menahem 
Department  of  Applied  Mathematics 
Weizman  Institute  of  Science 
Rehovot,  ISRAEL  951729 


Dr.  Robert  North 

Geophysics  Division 

Geological  Survey  of  Canada 

1  Observatory  Crescent 

Ottawa,  Ontario,  CANADA  K1A  0Y3 

Dr.  Frode  Ringdal 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjeller,  NORWAY 


Dr.  Jorg  Schlittenhardt 

Federal  Institute  for  Geosciences  &  Nat'l  Res. 
Postfach  510153 

D-3000  Hannover  51,  FEDERAL  REPUBLIC  OF 
GERMANY 


Prof.  Daniel  Walker 
University  of  Hawaii 
Institute  of  Geophysics 
Honolulu,  HI  96822 


Ms.  Eva  Johannisson 
Senior  Research  Officer 
National  Defense  Research  Inst. 
P.O.  Box  27322 
S- 102  54  Stockholm,  SWEDEN 
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